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ABSTRACT 


This  report  discusses  the  autonomous  orbit  control  system  launched  aboard 
the  synchronous  communications  satellite  LES-6  in  late  1968.  The  topics 
presented  include  the  motivation  for  automatic  orbit  control,  the  measure¬ 
ment  techniques,  observability  and  controllability,  the  design  of  a  variable 
limit-cycle  system,  automatic  station  acquisition,  simulation,  and  testing. 
The  history  of  the  development  of  this  control  system  is  emphasized,  and 
the  early  flight  results  are  discussed. 
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AUTOMATIC  ORBIT  CONTROL 
OF  THE  LINCOLN  EXPERIMENTAL  SATELLITE  LES-6 


I.  INTRODUCTION 

For  several  years,  satellites  have  been  used  to  establish  radio  communication  links  between 
stations  located  on  the  surface  of  the  earth.  Radio  signals  transmitted  by  one  station  are  re¬ 
ceived  by  the  satellite,  amplified,  and  retransmitted  to  one  or  more  terminal  points  on  the 
ground.  Normally,  line-of-sight  transmission  is  required,  so  that  only  those  stations  which 
would  be  visible  to  a  hypothetical  observer  on  the  satellite  can  be  connected  by  this  type  of  radio 
link. 

A  major  requirement  for  any  communications  system  is  that  the  link  between  stations  be 
uninterrupted,  so  that  transmission  of  information  between  the  terminals  is  possible  at  all  times. 
This  need  has  led  to  the  introduction  of  "synchronous"  satellites  which  revolve  around  the  earth 
at  an  average  angular  rate  equal  to  that  at  which  the  earth  spins  around  its  polar  axis.  As  is 
well  known,  the  orbital  period  of  a  satellite  is  related  to  its  total  energy,  which  in  turn  is  a  func¬ 
tion  of  the  orbital  major  axis.  By  appropriately  choosing  the  latter,  the  satellite  period  can  be 
made  equal  to  one  sidereal  day.  Furthermore,  if  the  orbit  is  equatorial  and  circular,  the  sat¬ 
ellite  is  geostationary,  or  motionless,  relative  to  an  observer  on  the  earth.  It  follows  that  one 
way  of  avoiding  interruption  of  transmission  due  to  the  motion  of  the  satellite  with  respect  to  the 
ground  terminals  is  to  place  the  spacecraft  in  a  circular,  equatorial,  synchronous  orbit. 

Unfortunately,  this  ideal  situation  cannot  be  achieved  permanently.  A  synchronous  satellite 
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is  subject  to  forces  which  tend  to  displace  it  from  the  geostationary  position:  a  stable  syn¬ 

chronous  orbit  exists  at  only  two  earth  longitudes,  about  107°  West  (over  the  Pacific  Ocean)  and 
73°  East  (over  the  Indian  Ocean),  In  general,  the  satellite  has  a  long-term  pendular  motion  about 
the  initial  synchronous  position,  which  must  be  opposed  by  a  thruster  on  the  satellite.  On  the 
first  synchronous  satellites,  the  thruster  was  commanded  by  a  signal  from  the  ground,  after  a 
difference  between  the  actual  and  the  desired  satellite  positions  had  been  detected  by  tracking, 
until  the  desired  drift  velocity  and  position  were  obtained, 

A  completely  automatic  and  self-contained  technique  has  been  realized  for  the  first  time  on 
the  Lincoln  Experimental  Satellite  LES-6,  a  communications  satellite  that  was  launched  success¬ 
fully  on  26  September  1968.  In  this  system,  the  angular  position  error  is  measured  on  the  satel¬ 
lite  itself  and  thrusters  are  fired  according  to  a  correction  sequence  determined  by  a  digital  con¬ 
troller  on  board.  This  approach  (a)  eliminates  tracking  and  commanding,  and  (b)  is  potentially 
more  accurate  than  the  ground -command  method.  The  satellite  can  also  be  made  to  acquire  new 
longitudinal  positions  by  changing  an  internal  reference,  so  that  it  can  follow  new  stations  in  a 
preprogrammed  sequence  or  upon  receipt  of  simple  instructions  from  the  ground. 

The  use  of  an  Automatic  Orbit  Control  system  significantly  reduces  the  need  for  ground 
tracking.  We  believe  that  this  is  a  primary  economic  incentive  to  the  use  of  these  techniques  in 
the  future.  To  appreciate  the  fact,  it  is  instructive  to  consider  our  own  experience  in  determin¬ 
ing  the  orbit  of  LES-6,  using  a  30-foot  UHF  antenna  in  Lexington,  Massachusetts.  The  data  re¬ 
corded  by  the  ground  station  are  satellite  range,  obtained  from  the  round-trip  time  of  a  signal 
sent  through  the  satellite  repeater,  and  elevation  and  azimuth  angles,  derived  directly  (and 
rather  inaccurately)  by  pointing  the  antenna  for  maximum  received  signal.  These  measurements, 
which  last  several  minutes,  are  repeated  hourly  for  30  consecutive  hours  and  are  then  processed 
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by  an  Orbit  Determination  Program,  which  performs  a  least-square  fit  on  the  data  to  yield  one 
value  of  satellite  position  and  drift.  This  orbit  determination  is  thus  a  non-trivial  burden,  which 
of  course  becomes  more  significant  as  the  number  of  satellites  to  be  tracked  increases,  or  when 
the  ground  station  must  be  used  for  purposes  other  than  satellite  tracking.  Ground  station  sat¬ 
uration  is,  therefore,  one  of  the  problems  that  Automatic  Orbit  Control,  done  directly  from  the 
satellite,  will  alleviate. 

The  second  advantage  of  an  Automatic  Orbit  Control  system,  i.e.,  that  of  accuracy,  is  still 
to  be  proven.  In  fact,  one  major  aim  in  launching  the  LES-6  system  has  been  to  collect  infor¬ 
mation  that  will  enable  us  to  design  high-accuracy  second  generation  systems.  Our  goal  on 
LES-6  was  a  rather  modest  ±2  degrees  in  steady-state  longitude  accuracy  this  is  amply  suffi¬ 
cient  for  earth-coverage  microwave  systems,  and  is  indeed  very  satisfactory,  as  is  known  to 
those  who  have  experience  in  synchronous  satellite  orbit  control.  The  longitudinal  accuracy 
which  can  be  achieved  with  ground  tracking  techniques  depends  on  (a)  the  orbit  determination 
accuracy,  and  (b)  the  time  interval  between  corrections.  The  Syncom  satellites,  for  example, 
have  achieved  accuracies  of  2  degrees,  with  corrections  about  every  30  days.  If  the  orbit  is 
measured  and  controlled  directly  from  the  satellite,  the  constraints  imposed  on  the  correction 
interval  by  ground  station  limitations  are  removed.  Greater  orbit  control  accuracy  is  important 
for  satellites  using  narrow-beam,  i.e.,  high-gain,  antennas.  A  second  incentive  for  accurate 
longitude  control  is  connected  to  the  use  of  many  synchronous  satellites  closely  spaced  in  longi¬ 
tude:  in  the  future  there  will  be  a  problem  of  allocation  of  "synchronous  orbit  space,"  especially 
in  the  neighborhood  of  some  preferred  longitudes.  In  this  report  we  discuss  some  preliminary 
statistical  results  obtained  from  sun  and  earth  sensors. 

The  work  that  is  presented  here  refers  to  near-circular  orbits,  and  the  flight  data  were  col¬ 
lected  from  the  synchronous  equatorial  satellite  LES-6.  However,  these  ideas,  and  those  which 
will  develop  from  our  own  and  others'  experience,  will  probably  be  applied  to  other  types  of  or¬ 
bits.  Earth-resources  satellites  and  large  unmanned  orbital  platforms  may  become  a  reality. 

As  the  number  of  active  satellites  increases  and  improvements  in  solar  power  systems  lengthen 
their  useful  lives,  on-board  control  of  the  orbit  parameters  will  become  increasingly  important. 
The  acceptance  of  automatic  techniques  depends,  of  course,  on  their  reliability.  At  the  time  of 
this  writing,  the  Automatic  Orbit  Control  system  of  LES-6  had  operated  properly  for  eleven 
months.  During  this  time,  the  Pulsed  Plasma  Thrusters  operated  for  about  3,500  hours.  We 
hope  to  be  able  to  demonstrate  that  the  system  can  control  the  orbit  automatically  for  several 
years . 

Automatic  orbit  control  has  not  been  attempted  until  now,  although  self-(  ontained  methods 
have  been  proposed  in  the  literature  by  several  authors.  Some  of  these  studies  include  applica¬ 
tions  of  filtering  theory,  especially  Kalman's,  to  this  problem.  The  explanation  of  the  delay  in 
realization  lies,  at  least  in  part,  in  the  complexity  of  the  orbit  determination  techniques  which 
had  been  proposed  until  recently.  Reliability,  weight,  and  power  limitations  have  made  all  the 
systems  proposed  impractical  for  today's  communications  satellites. 

It  is  interesting  to  review  briefly  the  development  of  self-contained  navigation  ideas,  which 

naturally  center  around  the  problem  of  determining  the  orbit  from  autonomous  data.  A  good  gen- 
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eral  survey  up  to  1965  is  made  by  Gansler,  who  considers  the  classic  distinction  between  pre¬ 
liminary  determination  and  updating  of  the  orbit ^  The  existing  literature  on  the  subject  of 

j  Gansler 's  evoluatian  of  the  characteristics  of  o  "typicol  "  sel  f-cantai  ned  navigotian  system,  excluding  the  pro¬ 
pulsion  equipment,  is  that  it  will  weigh  87  pounds  and  require  218  watts.  The  LES-6  system  weighs  4  pounds 
and  uses  120  mil  I  iwatts. 
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on-board  orbit  determination  falls  in  three  categories:  (a)  discussions  of  instruments  with  which 
to  observe  satellite  position  and  dynamic  coordinates,  (b)  elaborations  of  the  orbit  equations, 

(c)  applications  of  filtering  theory,  especially  Kalman's.  The  first  group  includes  the  many  types 
of  sensors,  star  trackers,  inertial  platforms,  and  dynamometers  which  have  been  proposed  but 
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will  not  be  discussed  here.  The  best  paper  in  the  second  group  is  by  Gersten  and  Schwarzbein, 
who  give  formulas  for  the  orbital  elements  as  a  function  of  various  combinations  of  measure¬ 
ments  of  satellite  altitude,  anomaly,  and  velocity  components.  The  difficulty  is  that  the  rela¬ 
tionships  are  transcendental  equations  which,  even  if  expressed  in  closed  form,  require  relatively 
sophisticated  computation  systems  to  obtain  reasonable  accuracy.  The  same  authors  examined 

Q 

some  additional  cases  in  which  intervals  of  time  as  well  as  dynamic  coordinates  are  measured; 

they  suggest  that  the  equation  of  motion  be  solved  repeatedly  for  different  values  of  orbital  major 
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axis  until  the  time  measurements  agree  with  the  computed  results.  McAllister  and  Wagner 

give  the  equations  of  motion  for  impulsive  perturbations  of  the  dynamic  coordinates  of  a  satellite 

in  a  circular  orbit  as  a  state  transition  matrix.  Several  authors  apply  Kalman's  statistical  fil- 
H_20  13 

tering  theory  to  the  determination  of  orbits.  Battin  also  solved  the  problem  of  selecting 
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the  best  star  to  use  as  a  reference  in  making  angular  measurements.  Frazier,  et  al.f  examined 
the  possibility  of  approximating  the  values  of  satellite  anomaly  and  the  subtense  of  the  principal 
body  with  polynomials  in  time  and  found  that,  in  general,  it  is  preferable  to  use  the  linearized 
form  of  the  known  dynamic  equations  to  smooth  the  data,  leading  to  more  complex  data  process¬ 
ing  requirements. 

In  196  5,  W.  E.  Morrow  of  Lincoln  Laboratory  proposed  that  a  simple  combination  of  sun 
and  earth  sensors  could  be  used  to  determine  the  satellite  "high -noon"  position  in  orbit;  a  stable 
oscillator  could  furnish  an  appropriate  reference  for  the  motion  in  sidereal  space  of  any  earth 
meridian,  and  the  correction  for  sun  transit  time  effects  could  be  performed  digitally  on  board. 
Shortly  thereafter,  the  author,  while  studying  the  problem  of  the  automatic  pointing  of  a  com¬ 
munication  antenna  from  a  near-circular  orbit  to  earth,  proposed  simple  methods  of  orbit  detcr- 
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mination  and  reconstruction  which  were  particularly  suited  for  implementation  on  a  satellite. 

The  present  LES-6  Automatic  Orbit  Control  system  has  evolved  from  these  ideas. 

Section  H  of  this  report  discusses  the  longitudinal  perturbations  of  synchronous  satellite 
orbits.  Section  HI  deals  with  the  orbit  determination  problem,  for  small  eccentricity  orbits; 
solutions  requiring  linear  combinations  of  times  of  flight  are  discussed.  Section  IV  examines 
the  observability  and  controllability  issues  from  the  point  of  view  of  modern  control  theory. 

Section  V  describes  the  finite-state  controller  which  was  actually  flown  on  LES-6  and  discusses 
the  automatic  station -acquisition  maneuver.  Section  VI  deals  briefly  with  the  hardware  imple¬ 
mentation,  the  simulation  studies,  the  choice  of  the  integration  routines  and  step  size,  the  ran¬ 
dom  and  deterministic  error  simulation,  and  the  testing  of  the  flight  system.  Section  VII  dis¬ 
cusses  the  dynamic  behavior  in  orbit  and  several  statistical  hypotheses  on  sensor  errors. 

II.  EFFECTS  OF  EQUATORIAL  ELLIPTICITY  ON  SYNCHRONOUS  ORBITS 

The  influence  of  the  asymmetries  of  the  earth's  gravitational  potential  on  the  motion  of  sat¬ 
ellites  has  been  studied  by  several  authors.  For  equatorial  synchronous  satellites,  the  most 
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important  effect  is  that  due  to  the  "equatorial  eccentricity"'  which  is  of  the  order  of  10  .  We 

f  One  should  consider  thot  this  is  simply  o  convenient  woy  of  interpreting  the  measured  gravitational  potential  in 
the  equatorial  plane.  In  foct,  far  an  eccentricity  af  10”^,  the  difference  between  mojar  ond  minor  semiaxes  is 
only  about  380  feet. 
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Fig.  1.  Reference  system  for  elliptic  earth. 


collect  here  some  of  the  important  results,  fol¬ 
lowing  mainly  Ref.  1.  We  further  develop  the 
"weak  coupling"  approximation  and  show  it  to 
be  valid  for  the  satellite  motion  sampled  once 
per  day. 

Stable  synchronous  orbits  exist  only  at  the 
two  longitudes  which  correspond  to  the  minor 
axis  of  the  equatorial  ellipse.  The  longitudes 
along  the  major  axis  are  unstable  equilibrium 
points.  At  all  other  longitudes,  the  satellite  is 
subject  to  a  drift  acceleration  and  its  motion 
exhibits  a  long-term  pendular  motion,  depend¬ 
ent  on  the  initial  conditions,  about  a  stable  point 
According  to  recent  determinations,  the  minor 
axis  is  located  at  about  73°  East.  Let  r,A  be 
defined  as  in  Tig.  1,  and  let  O  be  the  satellite 
colatitude.  Further,  R  =  mean  equatorial  ra¬ 


dius  of  the  earth,  =  coefficient  of  the  second  zonal  harmonic  of  the  earth  (which  is  related  to 
the  polar  oblateness),  and  y  ~  ellipticity  of  the  equator. 

The  kinetic  energy  in  a  coordinate  frame  which  rotates  rigidly  with  the  earth  at  angular 
speed  cl>  is 


_  1  -2  ,  1  2-2  ±  1  2  .  2^  ,•  v2 
T  -  ^  r  +  2re  +  2  r  sin  e  ^  +  w) 


and  the  potential  energy  can  be  written 
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We  wish  to  find  stationary  solutions  r  -  r  ,  A  A  ,  G  Q  .  From  Lagrange’s  equations, 

&  4^(t-v)  =  4(T~V)  (,["3) 

(where  q^  is  any  of  the  dynamic  coordinates  r.  A,  G)  follows,  for  the  stationary  case, 

sr-  (T  -  V)  =  0  , 

3qk 


and  hence  the  equations 
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The  last  equation  has  the  solutions  0^  -  0  and  7r/2,  of  which  the  last,  corresponding  to  equa¬ 
torial  orbits,  is  acceptable.  (This  can  be  seen  on  physical  grounds,  or  by  introducing  the  solu¬ 
tion  -  0  in  the  previous  equations  and  checking  that  no  solution  exists  for  rQ.)  Then,  from 
( II- 5) ,  stationary  points  are  found  to  exist  at  longitudes  0,  7r/2,  tt,  (3/2)  ir  (measured  from  the  ma¬ 
jor  axis  of  the  equatorial  ellipse).  The  values  of  rQ  which  correspond  to  these  solutions  are 
found  from 
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(the  positive  sign  corresponds  to  the  points  on  the  ellipse  major  axis).  To  discuss  the  stability 
of  these  solutions,  we  write  Lagrange's  equations  and  linearize  about  the  stationary  points. 

The  following  three  equations  result  easily: 
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( II -10) 


in  which  the  subscript  6  indicates  that  the  variables  are  variations  around  the  stationary  values. 

The  third  equation  is  decoupled  from  the  first  two.  The  motion  in  0  is  harmonic  with  fre¬ 
quency  oo/Zn,  i.e.,  with  daily  period,  given  the  orders  of  magnitude  of  the  quantities  involved. 

To  discuss  the  stability  of  the  other  two  coupled  equations,  we  find  the  eigenvalues  of  the 
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system,  in  which  we  can  make  approximations  remembering  that  p/c o  =*  r^  (this  can  be  obtained 
by  equating  the  centrifugal  and  gravitational  accelerations  in  synchronous  orbit)  and  that  J  and 
y  «  1.  Then  the  system  reduces  to 


*6  ~  3a,2r6  ~  2V*6  =  0  -  (II"11) 

I -f6  +  V2l'(r)2“\-°  •  <n-12) 

o  \  o/ 

Using  operator  notation,  the  equations  for  the  system  arc 


(p2  —  3 cd2)  z  ^  —  2rQCL)pz2  =  0  ,  (11-13) 

pz^  +  (p2  ±  ca>2)  z^  77  0  ,  (11-14) 

o 

2 

where  c  Zy  (R/rQ)  ,  from  which  the  eigenvalues  are  the  solutions  of  the  characteristic  equation 
p4  +  c2p2  *  3co;4  -  0  (11-15) 


5 


or 


P 


2 


-c2* 


\  1  ±  12c 


2 


-t3co;  =  ±6y 


(11-16) 


The  first  solutions  (±jcu)  are  oscillatory.  The  second  set  implies  instability  when  positive 
sign  applies  (major  axis)  and  stability  in  the  other  case  (minor  axis).  The  solutions  then  include 
a  daily  component  (at  frequency  co )  and  a  slow  oscillation  with  period 


T 


6.63 

Jg  x  10~S 


=*  850  days 


By  solving  the  system  with  appropriate  initial  conditions,  Blitzer,  et  al.,  find  that  the  trajectory 
in  the  rotating  coordinate  frame  is  a  very  elongated  ellipse  about  the  stable  longitude.  The  same 
study  shows  computer  results  for  the  complete  (nonlinear)  equations  of  motion.  Since  that  anal¬ 
ysis  was  made,  the  observations  of  synchronous  satellites  have  confirmed  its  validity.^ 

It  is  interesting  to  discuss  a  further  approximation  to  the  system  [  ( 1 1  —  8 ),  (II-9),  (11-10)], 
which  we  shall  call  ’’weak  coupling”  and  which  is  very  important  for  the  applications.  If  we  as¬ 
sume  r.  small,  we  can  write 
o 


3cor*  +  2r  A  "  -  0 
o  o  o 

- 

(II-17) 

—  rt  +  A  *  ±  Zy 
rod 
o 

/  R  \2  2.  *  A 

(_)  ,Vo  • 

(11-18) 

From  the  first  equation  it  is 

rJ  __  _  z  h 

r  3  a) 

(11-19) 

o 


This  equation  can  be  obtained  from  first  principles,  namely,  from  Kepler's  third  law.  It  relates 
small  changes  in  period  to  small  changes  in  semimajor  axis.  Taking  the  time  derivative  of 
(11-17)  and  substituting  in  (11-18), 


Gy 


(11-20) 


Thus,  we  have  obtained  a  second-order  equation  in  the  variable  A^  which  is  decoupled  from  the 
equation  in  r^  and  which  has  the  correct  long-term  motion  with  frequencies  ±  n/ Gy  (R/r  )  co  in 
the  stable  case. 

One  is  then  led  to  consider  what  happens  if  A^  is  sampled  once  a  day,  i.c.,  if  the  daily  com¬ 
ponent  in  the  solutions  is  desregarded.  In  order  to  check  that  indeed  the  variables  A^,  r^  of  the 
weak- coupling  approximation  are  the  values  of  A^,  r^  sampled  once  per  day,  wc  can  proceed  as 
follows. 

The  general  solution  for  the  linearized  system  [Eqs.  (11-11)  and  (11-12)]  in  the  stable  case  is 
obtained  easily 


f  It  must  be  noted  that  luni-solar  perturbations  alter  this  picture  by  effectively  changing  the  synchronous  semi¬ 
major  axis.  Closed  arbits  na  longer  exist  even  at  the  "stable"  paints. 
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■jWt]  -j  -4  ^  R  [B^™1  - 

n/6  1 


B2e 


-  j  «  a)  t , 


(11-22) 


where 


o'  -  \r6y  —  «  1 
o 

Let  us  consider  only  the  values  at  times  t  =  nTg,  where  Tg  2n/oo,  n  ~  0,  1,  2,.  .  .  .  In  particular, 
we  can  select  the  initial  value  (that  is  to  say,  the  phase  of  the  sampling)  so  that 

A,  -  A  -  0  .  (II-23) 

1  2 

Consider  then  the  values  of  r^,  r^,  and  A^  at  the  sampling  times,  which  we  indicate  with  a  star: 


r*  -j  J.  ^R[B.eJ"2,rn-n-e-J“2,rn|  ,  (11-24) 

6  n/6  1  Z 

r*  _L  Vy  R|«w)2  (bJ0,2111  -  B,e'j®2™|  ,  (11-25) 

6  '  \[F  1 

>*6-.  jau  (B1ej"27rn  -  B2e-i*2,rn]  .  (11-26) 

Remembering  the  definition  of  o\  we  obtain  that  oj^r'^  and  ojA'^  are  of  order  sfy,  while  r^  is  of 
order  [\ly)^  and  therefore  negligible  in  (II-8).  Thus,  we  can  conclude  that  if  one  uses  as  var¬ 
iables  the  radius  and  the  geocentric  longitude  observed  once  per*  day,  one  can  study  the  satellite 
drift  in  a  neighbor  hood  of  the  minor  axis  from  the  simple  decoupled  Eqs.  (11-19)  and  (11-20). 

The  previous  equations  and  conclusions  apply  to  motion  in  a  neighborhood  of  the  stationary 
points.  It  is  important  to  establish  the  weak-coupling  approximation  for*  all  values  of  longitude  A. 
Let  us  linearize  about  the  values  r  and  all  derivatives  equal  zero,  treating  A  as  a  finite  quan¬ 
tity.  Then  Kq.  (II -8)  remains  unchanged,  while  (11-9)  becomes 


2oj 

r 
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4  A  4 


ypR 


>2 


sin  2A 


0 


(11-27) 


an  equation  which  is  nonlinear-  in  A.  Let  us  consider  linearized  solutions  to  the  previous  equa¬ 
tion.  For  A  ~  A  we  can  write 
o 


Vj.  -  3a;2rx  -  2r  ojA ,  0 

6  6  o  6 

^  +XA  =-  IE”!  sin 2A 

v  6  6  5  o 

o  r 

o 

The  characteristic  equation  is  found  easily, 

4  .  2  2  n 

p  +  co  p  0  , 

and  the  general  solution  is 


(11-28) 

(11-29) 


(11-30) 
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X6(t)  =  Al6' 


(11-31) 


+  A-.e_ja;t  +  B.  +  B_t  +  4  ^5-  sin  2A  '•  t2 


2  2  5 

r 

o 


r.  (t)  =  —  A.  -y-  je*^*  +  A?  ~~  je”‘®a;t  —  2  —  — — sin 
6  1  2  J  2  2  J  c o  5 


?  r 

2A  •  t  -  4  —  B7 
o  3  cj  2 


Sampling  at  t  -  nTg  <2n/w)  n, 

2 

A*(nT)  =  (A  +  A2  +  B1 )  +  B2nT  +  sin2Af)  •  (nT)2 

r 

o 

Thus  the  longitude  sampled  once  per  day  satisfies  the  differential  equation 

A*  -  sin  2A  =  0 

6  5  o 

r 

o 


(11-32) 


(11-33) 


(11-34) 


We  can  write  the  general  solution  to  (11-27)  by  letting  the  coefficients  of  the  linearized  solu¬ 
tion,  Eq.  (11-31),  be  functions  of  time: 

A*(t)  =  DQ  +  U^t)  •  t  +  D2(t)  •  t2  ,  (H-35) 

where  the  coefficients  D^(t),  D^ft)  are  determined  by  choosing  different  values  for  X Then  we 
ean  conclude  that  the  weak-eoupling  approximation  is  valid  for  the  nonlinear  equation 

A*(t)-  sin  2A*  =  0  .  (11-36) 

r 

o 

This  is  the  well-known  “pendulum  equation"  of  elementary  dynamics.  For  small  excursions 
of  A  ,  the  linearized  approximation,  leading  to  parabolic  motion,  Eq.  (11-33),  can  be  used.  As 
we  shall  see,  the  pendulum  equation  is  very  useful  in  practice,  since  the  motion  of  the  satellite 
is  thus  reduced  to  a  one-dimensional  problem.  The  radius  vector  is  still  given  by  (11-19).  It  ean 
further  be  shown  that  considering  the  higher-order  tesseral  harmonics  does  not  essentially 
modify  these  conclusions. 


III.  SIMPLE  METHODS  TO  DETERMINE  NEAR-CIRCULAR  ORBITS 
A.  General  Principles 

Since  elliptic  orbits  ean  be  described  by  any  set  of  sufficient  parameters,  we  wish  to  use 
parameters  which  are  most  easily  measured  with  on-board  instrumentation.  The  orbit  can  then 
be  controlled  by  maintaining  these  parameters  at  a  reference  value:  the  computation  problem, 
and  the  measurement  and  control  problems  are  treated  as  a  whole.  The  desired  control  accu¬ 
racy  determines  the  precision  necessary  in  the  measurement  and  in  the  computations,  and  this 
in  turn  affects  the  choice  of  the  technique  of  orbit  determination.  The  methods  summarized  here 
only  require  linear  combinations  of  measured  times  of  flight.  The  results  arc  obtained  from 
well-known  expansions  of  the  law  of  motion  in  terms  of  the  orbital  eccentricity  e  to  terms  of  the 
third  order.  Several  properties  of  the  seeond-order  solution,  which  are  of  great  practical  inter¬ 
est,  are  discussed.  These  techniques  can  also  be  applied  to  the  determination  of  the  center  of 

2 1 

the  earth  from  the  satellite  and  to  circularize  an  orbit  automatically,  as  seen  in  the  following. 
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Sinee,  as  we  have  seen,  the  out-of-plane  perturbations  of  the  motion  of  the  satellite^  are  de¬ 
coupled  from  the  motion  in  the  orbit  plane,  we  study  plane  motion  only. 

The  instrumentation  needed  for  the  measurements  eonsists  of  an  earth  sensor,  a  sun  or  star 

sensor,  and  a  clock.  The  earth  sensor  must  determine  the  direction  of  the  center  of  the  earth, 

2  3-25 

or  loeal  vertieal,  and  ean  be  an  infrared  or  optieal  sensor.  If  an  optieal  sensor  is  used, 

electronic  processing  of  its  output  is  necessary  to  correct  for  the  variation  of  apparent  width  of 
the  earth  as  the  earth- satellite- sun  angle  ehanges.  The  sensor  is  relatively  simple  if  the  orbit 
is  nearly  eireular,  so  tnat  the  earth's  angular  subtense  varies  little  in  time,  and  if  the  semimajor 
axis  equals  a  few  earth  radii.  The  determination  of  the  loeal  vertical  when  these  two  conditions 
are  not  satisfied  requires  more  sophisticated  sensors.  The  difficulty  of  this  measurement  per¬ 
mits  one  to  state  that  the  precision  of  orbit  determination  from  the  satellite  is  at  present  essen¬ 
tially  determined  by  the  aeeuraey  with  whieh  the  loeal  vertieal  can  be  established.  Several  types 
of  sun  and  star  sensors  have  been  used  or  proposed,  the  latter  being  potentially  more  aeeurate, 
but  more  eomplex  than  the  former.  Star  sensors  are  generally  used  in  a  continuous  tracking 
mode,  rather  than  discontinuously  on  spinning  satellites.  The  advantages  of  using  star  trackers 
are  that  the  angular  subtense  of  stars  is  much  smaller  than  the  sun's  and  that  for  all  praetieal 
purposes  they  can  be  considered  a  true  inertial  reference.  The  clock,  a  very  aeeurate  and  stable 
oscillator,  in  general  uses  a  temperature-eontrolled  quartz  crystal;  long-term  stability  of  one 
part  in  10  million  ean  be  obtained  on  today's  satellites. 

The  planes  of  sight  of  both  the  earth  and  sun  sensors  are  typieally  perpendicular  to  the  or¬ 
bital  plane,  eontain  the  spin  axis  of  the  satellite,  and  form  an  angle  6.  As  the  satellite  spins, 
the  sun  and  earth  sensors  deliver  two  trains 
of  pulses  which  in  general  do  not  occur  at  the 
same  time.  However,  for  any  given  orbit,  an¬ 
gular  separation  of  the  sensors,  and  position  of 
the  sun,  there  is  one  point  P  in  orbit  at  whieh 
the  pulse  trains  from  the  two  sensors  are  in 
phase  (Fig.  2).  If  another  sun  sensor  is  used 
together  with  the  earth  sensor,  at  an  angular 
separation  <5\  a  seeond  point  P'  ean  be  estab¬ 
lished  in  orbit.  Alternatively,  one  eould  use 
two  earth  sensors  and  one  sun  sensor.  Thus 
each  combination  of  three  sensors  is  an  orbital 
angle  finder,  which  together  with  a  elock  can 
measure  the  time  in  whieh  the  satellite  sweeps 

A 

the  known  angle  PEP'.  If  the  sun,  rather  than 
a  star,  is  used  as  a  reference,  its  angular  mo¬ 
tion  relative  to  the  orbit  must  be  taken  into  aeeount.  In  the  case  of  an  elliptic  orbit,  the  argu¬ 
ment  (or  angular  coordinate)  of  the  points  P  and  P',  relative  to  perigee,  changes  day  by  day 
beeause  the  projection  of  the  sun  line  on  the  orbital  plane  moves  with  yearly  period  and  beeause 
of  perigee  preeession.  The  orbital  angle  finder  and  the  cloek,  in  measuring  the  time  of  flight  of 
the  satellite  between  points  in  orbit,  essentially  measure  the  average  angular  velocity  of  the  sat¬ 
ellite  as  it  sweeps  angles  of  known  width. 

t  It  is  well  known  that  the  orbital  plane  of  an  earth  satellite  is  perturbed  mainly  by  luni-salar  effects  and  by  the 
spherical  harmonics  In  the  geapatential .  The  magnitude  af  these  perturbations  depends  an  the  arbit  considered. 


~]~~3-63-a044T] 


Fig.  2.  Measurement  af  times  af  flight. 
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(111-1) 


Expanding  the  law  of  motion  in  powers  of  the  eccentricity  e, 

wt  =  e  -eo  —  2c  [sin  (©)  —  sin  (0Q)]  +  |  c2  (sin(20)  -sin(20o)] 

where  0  =  satellite  argument  measured  from  perigee,  =  orbital  radian  frequency  0=0^  for 
t  =  0. 


Fig.  3.  Measurement  points 
for  orbit  reconstruction. 


Consider  the  times  in  which  the  satellite  sweeps  three  consecutive  90°  sectors  in  orbit,  de¬ 
fined  by  points  P^,  P^,  P^,  and  P^  (Fig.  3),  and  let  0q  be  the  argument  at  P^.  The  satellite 

therefore  moves  from  P,  to  P.  in  a  time  interval  t,  such  that 

12  1 

wt^  =  j  —  2e  [sin  (0Q  +  j  )  —  sin  (0q )]  +  e2  [sin  (20q  +  7r)  -  sin  ( Z Oq ) ]  .  (Ill- 2) 

Writing  similarly  the  time  intervals  t^,  t^  in  which  the  satellite  travels  from  P^  to  P^  and  to 
P^,  respectively,  and  taking  linear  combinations  of  the  equations,  one  finds  the  (first-order) 
parameters: 

e  sin(Qo)  =  |  [w(t1  +  t2)  -  tt]  ,  (111-3) 

e  cos(Oo)  =  |  [tt-  o)(t2  +  t3)]  .  (Ill—  4 ) 

The  quantities  o>,  e  sin(0Q),  e  cos  (0  )  can  be  used  as  sufficient  parameters,  from  which 
e2  cos  (20  )  can  be  determined.  However,  since  this  computation  is  not  easily  implemented, 
e2  sin(20Q)  are  considered  additional  (second-order)  parameters  for  the  equation  of  motion;  this 
is  legitimate,  since  enough  independent  measurements  of  times  of  flight  can  be  made. 

To  obtain  the  second- order  parameters  independently,  the  time  t*  in  which  the  satellite 
travels  at  45°  sector  from  P^  to  P'  (Fig.  3)  is  measured,  being 

wt*  =  |  -  2e  lsin(0o  +  ^)- sin(eo)]-  |  e2  [cos(20o)-sin(20o)]  .  (HI- 5 ) 

Using  this  result  and  the  previous  equations,  after  some  simplifications,  the  reconstruction  of 
the  law  of  motion  in  terms  of  these  measurements  is 
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e  -e 


- 1 


2  I  "2  ~  ^2  +  l3M  sin  (^’l)  +  \  +  t2)  -  •  [cos  (cot )  —  1] 


5 

+  3 


3  Z~  1,  +  j  1  +  1  \ZZ  t  t* 

4  1  2  2  4  3 


*  sin(2o;t) 


+  ^  I  \  -  (t1  +  t3 )]  [cos  (2wt)  -  1 1 


(1II-6) 


to  third  order  in  €.  This  expression  involves  only  weighed  sums  and  differences  of  the  times 
of  flight,  as  coefficients  of  the  harmonic  terms. 

Some  simple  properties  of  the  second-order  reconstruction  of  the  orbit  are  of  practical  in¬ 
terest  because  satellites  are  often  placed  in  circular  orbits  around  the  principal  body.  Return¬ 
ing  to  (It  1-2)  and  dropping  the  terms  in  t  ^ , 


<  s  i  n  (0^ ) 


1_ 

4 


a*(t2  t3) 


,  cos(eo)  T^(t1  t2) 


(1 1 1-7 ) 


which  are  more  compact  forms  of  (HI-3)  and  (HI-4)  for  the  second-order  case.  These  expres¬ 
sions  can  be  used  to  rewrite  (III-6)  for  this  case.  The  orbit  reconst  ruction  is  now  particularly 
simple,  being 


e  -e 


- 1 


(t1  t2)  sin(cot)  + 


(t 3  ~  t2 )  [  1  -  cos  (cot)] 


(111-8) 


An  interesting  second-order  property  is  that  the  times  t',  tM  in  which  the  satellite  sweeps 
any  two  angles,  both  of  width  A0  and  opposed  with  respect  to  the  focus  are,  to  second  order  in  t, 

wt*  =  A0  -  Zi  [sin(0Q  +  AQ)  -  sin(0  )]  ,  (HI-9) 

ajtM  A0  +  2(  [sin(0Q  +  A0)  -sin(0Q)]  ,  (111-10) 

and,  therefore, 

T  7r  (t‘  +  tM)/AO  (111-11) 

Th i s  equation  provides  alternate  ways  of  measuring  the  period  for  quasi-circular  orbits  and  is 
particularly  simple  if  A0  is  an  integer  fraction  of  a  revolution.  To  sec  ond  order  in  c,  clearly 
it  is  also 


u  -  J  (u)'  4  a'")  ,  (H1--I2) 

where  co '  and  a? 1  ’  are  instantaneous  angular  velocities  at  two  points  in  orbit  opposed  with  respect 
to  the  focus. 

As  has  been  seen,  the  initial  condition  on  the  orbit  is  established  by  observing  the  time  of 
coincidence  of  the  outputs  of  a  sun  or  star  sensor  and  an  earth  sensor.  If  the  sun  is  used  as  a 
reference,  the  on-board  clock  must  be  a  solar  clock,  i.e.,  it  must  keep  into  account  the  valuation 
of  the  transit  time  of  the  sun  over  a  meridian  of  the  satellite.  If  the  orbital  reconstruction  is 
done  on  a  day-by-day  basis,  it  is  sufficient  for  some  applications  to  consider  the  mean  sun  mo¬ 
tion,  thus,  one  can  use  a  constant-speed  clock  and  correct  for  the  mean  sun  motion  either  in  car¬ 
rying  out  the  computation  or  by  altering  the  angular  distance  of  the  sensors  on  the  satellite  by  an 
appropriate  amount.  If,  on  the  other  hand,  the  initial  conditions  in  each  orbit  must  be  accurately 
known  for  several  years  (e.g.,  for  stationkeeping),  it  is  necessary  to  use  a  true  solar  clock, 
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which  can  be  implemented  electronically.  When  the  satellite  spins  about  an  axis  which  is  not 
parallel  to  the  earth's,  the  satellite's  true  solar  time  differs  from  the  earth's  (see  Appendix  A). 

In  determining  the  initial  conditions  on  elliptic  orbits  for  many  consecutive  days,  one  must 
consider  the  motion  of  the  line  apsides,  which  causes  the  point  of  coincidence  between  sensor 
outputs  to  move  day  by  day  along  the  orbit.  Furthermore,  if  the  sun  is  used  as  a  reference,  the 
point  of  coincidence  moves  with  the  projection  of  the  sun  line  along  the  orbit.  Thus,  if  the  time 
coincidence  occurs  at  point  P1  on  the  first  orbit,  it  occurs  at  P2,  P3,  .  .  .  Pn  on  successive  orbits 
Then,  the  n-th  coincidence  will  occur  at  a  time 

t  nT  +  At  (<p  ,  0  ) 
n  vvn  o' 

liven  though  a  true  solar  time  is  used,  the  fact  that  the  angle  is  traveled  on  an  elliptic  orbit 
is  cause  of  a  measurement  error.  The  error  varies  as 

6  2(  sin  (</>n)  ( t  ad) 

to  second  order  in  <  . 

If  the  measurement  is  continued  for  a  year  or  more,  so  that  tp  varies  from  0  to  2i r,  the 
peak-to-peak  error  is  6^  ^  -  228  <  (deg). 

A  simple  way  of  eliminating  this  measurement  error  to  second  order  in  e  is  then  to  repeat 
the  determination  of  the  initial  condition  twice  per  orbit  at  two  points  180°  apart  and  to  take  the 
arithmetic  average  of  the  two  determinations  j  This  technique  should  be  employed  also  when  a 
star  tracker  is  used  in  the  orbital  goniometer  to  eliminate  the  effect  of  apsidal  line  motion. 


B.  Application  to  Automatic  Orbit  Circularization 

As  an  application  of  these  ideas,  we  shall  describe  a  system  that  would  automatically  reduce 
e  $  The  controller  is  of  the  null -seeking  type,  requires  little  computation  on  the  satellite,  and 


Fig.  4.  Automatic  circularization 
technique. 


does  not  interact  with  the  control  of  the  satellite  daily  longitude.  The  orbital  position  is  meas¬ 
ured  at  points  P^  P^,  P^,  P^,  90°  apart  (Fig.  4),  with  a  system  of  angle-finders  of  the  type 
which  has  been  described  previously.  Pet  t^,  t^,  t^,  and  be  the  times  in  which  the  satellite 

t  U.S.  Patent  Application  643493,  filed  5  June  1967. 

$  This  system  has  not  been  implemented  an  LES-6. 
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moves  along  four  consecutive  90°  sectors.  Instead  of  computing  the  argument  of  apogee,  the 
quadrant  in  which  apogee  is  located  with  respect  to  the  measured  positions  can  be  determined  by 
noting  that  when  tj  =  max  t.  (i  1, 2,  3,  4),  apogee  is  between  P^  and  +  with  the  convention 
that  h  +  1  1  for  h  -  4,  In  the  boundary  case,  t  ^  >  t  ^  t^,  apogee  is  at  P^,  and  similar  relations 

hold  when  apogee  is  at  P^,  P^,  and  Thus,  the  position  of  apogee  is  determined  qualitatively 

without  computations,  using  only  comparisons  of  times  of  flight. 

The  longitude  (stationkeeping)  control  thrust  is  applied  symmetrically  about  the  orbital  focus 
so  that  the  eeeentricity  is  not  modified.!  For  eonvenienee,  if  F  is  the  impulsive  stationkeeping 
eontrol  thrust  along  the  orbit  tangent,  4/4  lr  shall  be  applred  four  times  per  orbit,  at  points  P^ 
through  P  .  To  change  the  orbital  eccentricity  without  affeeting  the  satellite  period,  the  net 
tangential  force  applied  during  one  orbit  must  be  zero.  A  eontrol  aetion  +  F  is  applied,  at  each 
of  the  two  fix  points  which  encompass  apogee,  followed  by  an  action  — F  at  the  two  points  encom¬ 
passing  perigee.  This  system  of  forces  ±T  tends  to  circularize  the  orbit  and  does  not  modify 
2  3  C 

the  satellite  period.  In  the  boundary  eases,  the  control  aetion  is  modified  simply.  The  time 
compar  isons  are  repeated,  and  the  corrective  action  ±F*e  is  applied  until  t.  -  tj  (i,  h  1,  2,  3,  4), 
within  the  limits  of  accuracy  of  the  sensors.  F  can  be  designed  to  satisfy  fuel  constraints  or 
other  requirements.  If  the  stationkeeping  corrections  1  and  the  circularizing  eorreet ion  ±F 
are  superimposed,  the  system  of  forces  shown  in  Fig.  4  results. 

To  determine  the  precision  of  this  technique,  t lie  differences  between  times  of  flight  are 
written  to  second  order  in  € 


-t2|  -  4c  I  cos(eo)| 


«h2  -t3!  4c  I  sin  (Oq)  | 


(111-13) 


The  control  system  operates  to  reduce  the  measured  differences  to  zero,  or 

w<h  -V  +CTi  u 2  0  • 


w(t2  ~  t^)  4-  0^2  -  W3  0 


(111-14) 


where  cj.  is  t lie  angular  error  made  in  measuring  the  corresponding  90°  sector. 

If  g.  is  further  written  as 

i 

cr .  ct  ^  +  (7^  ,  i  0,1,2  , 

it 

i.e,,  as  the  sum  of  a  zero-mean  and  a  systematic  component,  it  follows  that  this  control  tech¬ 
nique  is  insensitive  to  systematic  angle  errors.  It  is  interesting  to  see  what  e  one  would  have 
from  a  constant  measurement  error  cr  ^  -  a  ^  a ^  cr.  Then  from  (111-14), 

cr  4i  I  cos  (G  )|  ,  a  4(  I  sin  (6  ) 

r 1  o  1  r  o 


(in  rad),  or 

c  r  =  a/ 8“* 


f  This  statement  is  proved  in  Sec.  IV. 
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IV.  OBSERVABILITY  AND  CONTROLLABILITY  OF  NEAR-CIRCULAR  ORBITS 


In  this  section,  the  dynamic  problem  of  perturbations  of  near-circular  orbits  is  discussed 

from  the  point  of  view  of  modern  control  theory.  The  problem  of  observability  and  controllability 

27-29 

has  been  studied  by  several  authors  independently.  R.  Brockett's  results  are  the  most  com¬ 

plete  and  arc  used  in  parts  A  and  B  of  this  section.  In  part  C,  we  discuss  the  meaning  of  the 
loss  of  controllability  which  occurs  when  the  orbit  equations  are  sampled  once  per  day.  Wc  show 
that  the  symmetric  thrusting  policy  which  was  discussed  in  Sec.  Ill  does  not  change  the  orbital 
eccentricity.  This,  together  with  a  constraint  which  is  shown  to  exist  between  the  uncontrolled 
variables,  permits  us  to  show  that  this  thrust  policy  leads  to  a  Liapunov -stable  system  (i.e.,  one 
in  which  the  uncontrolled  variables  cannot  exceed  a  bound,  which  is  proportional  to  the  initial 
eccentricity).  Thus  we  prove  the  feasibility  of  controlling  the  satellite  longitude  using  observa¬ 
tions  appropriately  made  once  per  orbital  period,  and  reconcile  the  classical  approach  presented 
in  Secs.  11  and  111  and  the  point  of  view  of  modern  control  theory. 


A.  The  Continuous  Case 

Let  us  consider  the  orbit  equations  in  the  case  in  which  a  small  acceleration  u  is  imparted 
to  the  satellite  by  on-board  thrusters.  If  we  call  u^,  u^,  u^  the  three  components  along  the  dy¬ 
namic  coordinates  r,  r  A,  r^O,  the  linearized  equations  of  motion  [( 11-8),  (11-9),  (11-10)]  become 


i'*  3 to  r*  -  Zv  _ 

6  6  o  6 


r  A .  +  2w  r . 
o  o  o  <5 


r  G.  +  r  w  O. 
o  <5  o  5 


(I  V -1 ) 
( 1 V  -2 ) 

(IV -3) 


where  we  have  considered  the  spherical  components  of  the  earth's  potential  only,  since  we  can 

superimpose  the  perturbations  as  long  as  the  linearized  equations  hold. 

Taking  as  variables  x.  r*  x„  r.,  x_  r  A-  x. 
h  1  6*  Z  6*  3  o  6’  4 

form  of  the  resulting  differential  equation  is 


r  A  .,  x  _ 
o  6  5 


r  Or ,  x, 
o  o  6 


r  O  the  mat  r  ix 
o  (V 


x  ( t )  Ax(t)  +  Bu 


(1V-4) 


wit  h 


1  0  0  0  0 


0  0  2tu  0  0 


0  0  1 


0  0 


Zlx)  0  0  0  0 


0  0  0  0  1 


0  0  0  0  -cr  0 


(1V-S) 
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and 


0  0  0 

A  0  0 

0  0  0 

0  1  0 

0  0  0 

0  0  1 


One  ran  solve  the  equation  by  Laplace 


v 

1 

X1 

2  2 

S  -f  GO 

V 

~2co 

X3 

.  2  ,  2 , 
s(s  1  CO  ) 

xs 

0 

- 

transformation  and  find  with  easy  steps 


2c o 

0 

) 

/  2  ‘  2" 
s(s  i  W  ) 

1  i 

2  .  2 
s  —  3o) 

2.2  ,  2v 

S  (s  +  (JO  ) 

0 

U2 

0 

1 

T 

2  t  2 

S  t  O? 

(IV  -6) 


(IV  -7) 


The  natural  r  esponse  of  Kq.  (1V-4)  with  initial  values  x(0)  is  given,  in  the  time  domain,  by 

(IV  -8) 


x(t)  e—  x(0) 


Kt 


wher*e  e—  is  a  n  X  n  matrix. 


e-  I, "Vis  A)-1  T  i  (Is  A)-1eslds  , 

“  Z7T.1  Jp  -  “* 

I  being  the  identity  matrix  and  T'  a  contour  enclosing  all  the  eigenvalues  of  (2s  AJ* 
straightforward  but  lengthy  calculations  one  obtains 


4  t  cos  u 

sinuit 

U) 

0 

2(1  cos  ui ) 

U 

0 

0 

i  u>  sin  wl 

cos  ut 

0 

2  sin u.  l 

0 

0 

f>(- ut  +  sin  cjI ) 

-2(1  cos^t) 

(  lut  4  4  Kinu,'l) 

0 

0 

At 

e— 

CO 

U? 

6u(  It  cos  wt) 

2  sin  tot 

0 

(  t  4  4  cos  to  t ) 

0 

0 

0 

0 

0 

0 

cos  tot 

(sin  ut ) 

Zj 

0 

0 

0 

0 

— oj  smut 

cos  ut 

(IV  -9) 


After  some 


(1\  -10) 


IS 


The  solution  of  the  foreed  equation 


is  given  by 


x  (t )  =  Ax(t)  +  Bu(t) 


x(t)  =  e-(t)x(0)  +  f  e-(t-t7)  Bu(q-)  da 
'  o 


(IV-11) 


(1V-12) 


Let  us  consider  now  the  observability  of  the  linearized,  near-eireular  equations  of  motion 
in  the  continuous  ease.  Since  the  last  two  equations  in  (IV-4)  are  decoupled  from  the  others,  we 
can  treat  them  separately. 

As  is  well  known,  we  must  determine  the  rank  of  the  matrix 

L  (£.  ca, . . . ,  cAn_1) 

whose  rows  are  the  rows  of  C,  CA,  ete.,  where  Cx  is  the  observed  vector.  If  we  observe  the 
radius  only, 


Ci  (1, 0,0,0) 

and 


1  0  0  0 

0  10  0 

3  a;2  0  0  Zu 

0  -cj2  0  0 


(IV-13) 


vSinee  rank  [L^]  3  and  the  state  vector  x  is  of  dimension  4,  the  system  is  not  observable  by 

measuring  radius  (or  its  derivative)  only. 

Now  consider  looking  at  the  G  dependent  variable  x^.  Here  we  have 

C2  (0,  0,  1,  0) 

and 


0  0  10 

0  0  0  1 

0  -Zoo  0  0 


(IV- 14) 


—  6c 0  0  —  Aoo^ 


This  matrix  is  of  rank  4  so  the  system  is  observable  from 


16 


For  the  last  equation  we  have 


L 


5 


1 

0 


0 

1 


(IV-  15) 


which  is  nonsingular  and  therefore  ensures  the  observability  of  and  x^. 

Let  us  now  consider  the  controllability  of  the  system  described  by  Eq.  (IV- 4),  which  is  de¬ 
cided  by  considering  the  matrix 


K  --  (B,  AB,  .  .  An_1  B) 

If  wc  control  from  u.  (radial  thrust  alone), 

1 


0 

1 

0 

r~ - 

r\] 

3 

1 

1 

0 

2 

—  00 

0 

0 

0 

—  Zoo 

0 

0 

-  2oo 

0 

2u3 

Since  the  rank  of  this  matrix  is  only  3,  control  from  u^  is  not  possible. 

In  the  ease  where  u.  is  zero  and  u^  is  the  control,  the  K  matrix  takes  the  form 
1  2  — 


0  0  2w  0 

0  Zw  0  2  0> 

0  \  0  -4u>2 


0  —  4oo 


(IV'-  17 ) 


This  matrix  is  of  rank  4  and  hence  the  system  is  controllable  from  u^  alone.  Thus  we  see  that  if 
only  one  thrust  component  is  to  be  used  it  should  be  tangential  to  the  orbit. 

For  the  last  two  equations  we  have 


0  1 

1  0 


(IV-  18) 


which  is  clearly  nonsingular,  lienee  u^  can  control  the  out-of-plane  motion. 


B.  The  Sampled  Case 

As  we  have  seen  in  See.  Ill,  simple  combinations  of  sensors  permit  us  to  establish  points,  or 
fixes,  in  the  orbit.  Thus  we  are  interested  in  examining  what  happens  to  the  orbit  equations  when 
they  are  sampled  onee,  twice,  and  four  times  per  orbit.  We  also  assume  that  the  control  u  is 
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applied  only  at  the  sampling  times.  Then  u(kT)  is  the  magnitude  of  the  impulse  applied  at  (kT). 
Thus  the  system  is  deseribed  by  the  difference  equation  [from  (IV- 12)], 


x(kT  +  T)  =  e-(T)x(kT)  +  e-(T)Bu(kT) 


( IV— 19) 


where  T  =  sampling  period  =  2  7r/k a;,  (k  -  1,2,4).  The  transition  matrices  for  the  various  eases 

A(t ) 

are  easily  obtained  from  the  expression  of  e — v  '  (IV- 10).  These  are: 
i.  For  T  =  2t r/o;: 


X1 

1 

0 

0 

0 

0 

0 

X1 

0 

0 

0 

X2 

0 

1 

0 

0 

0 

0 

x2 

1 

0 

0 

X3 

(Tk+T)  = 

—  1  2  7T 

0 

1 

—  67T 

CO 

0 

0 

X  3 

(kT)  + 

0 

—  67 r 

CO 

0 

X4 

0 

0 

0 

1 

0 

0 

X4 

0 

1 

0 

X5 

0 

0 

0 

0 

1 

0 

X5 

0 

0 

0 

x6 

0 

0 

0 

0 

0 

1 

x6 

0 

0 

1 

_ 

_ 

ii.  For  T  =  7t/oj: 


(kT) 


( I V- 2  0 ) 


X1 

7 

0 

0 

4 

CO 

0 

0 

X1 

0 

4_ 

0 ) 

0 

X2 

0 

-  1 

0 

0 

0 

0 

x2 

-  1 

0 

0 

X  3 

(kT+T)  -- 

—  6n 

-  4 

CO 

1 

—  3  7T 

a; 

0 

0 

X  3 

(kT)  + 

-  4 

CO 

—  37T 

CO 

0 

X4 

—  12a; 

0 

0 

-7 

0 

0 

x4 

0 

-  7 

0 

x5 

0 

0 

0 

0 

-  1 

0 

X5 

0 

0 

0 

x6 

0 

0 

0 

0 

0 

- 1 

x6 

0 

0 

-  1 

- 

_ 

(kT) 


( I V  -  2 1 ) 


iii.  For  T  ~  tt / 2a; : 


X  A 

4 

1 

0 

2 

0 

0 

x . 

1 

2_ 

0 

1 

a; 

CO 

1 

a; 

a; 

X2 

3a; 

0 

0 

2 

0 

0 

X2 

0 

2 

0 

X. 

—  3tt  +  6 

-  2 

1 

(—  3  7T  +  8  ) 

0 

0 

x 

-2 

( —  3  7T  +  8) 

0 

3 

(kT+T)  = 

CO 

2a; 

3 

x  ( k  T )  + 

a ; 

2a; 

X  4 

—  6a; 

-2 

0 

-  3 

0 

0 

X4 

-2 

—  3 

0 

X5 

0 

0 

0 

0 

0 

1 

CO 

X5 

0 

0 

t_ 

CO 

x6 

0 

0 

0 

0 

—  a ; 

0 

x6 

0 

0 

0 

L  J 

- 

bJ 

* 

J 

u(kT) 


( I V—  22 ) 
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The  z-transform  of  the  responses  can  be  obtained  from  these  transfer  matrices.  Alterna¬ 
tively,  one  can  start  from  the  Laplace  transfer  function  (IV-7),  and  find  its  z-transform  for  the 
appropriate  value  of  T. 

Using  a  tilda  to  indicate  z-transform  functions,  one  obtains  easily: 


i.  For  T  =  Ztt/oj: 


0  0 


+  67TZ 


(x){Z  —  1  ) 


(IV-  2  3 ) 


0 


u 


3 


ii.  For  T  7r/o;: 


0 

4z 

(jO  (Z2  —  1  ) 

0 

U1 

—  4z 

+  3ttz 

0 

U2 

y2 

o;(z2  —  1 ) 

qj  (z  —  1  )2 

0 

0 

0 

- , 

m 

■ 

iii.  For  T  =  tt/Zq): 


z/(a) 

2  z(z  +  1 ) 

0 

“i 

yi 

2 

z  +  1 

W  (z  -  1  )  (z2  +  1  ) 

—  2z  (z  +  1 ) 

z  [(8  —  3 7T )  z2  —  l6z  +  (8  —  3 7T )] 

0 

U2 

y2 

Ci)(z  —  1 )  (z2  +  1  ) 

2o;(z2  +  1  )  (z2  —  1  ) 

0 

0 

z/a; 

^3 

y  3 

(Z2  +  1) 

(IV-24) 


(1V-25) 


The  test  for  observability  and  controllability  is  of  the  same  form  as  for  the  continuous  case. 
However,  the  A_  and  B  matrices  must  be  replaced  by  the  matrices  in  Eqs.  (IV-20),  (IV-21), 
(IV-22)  for  the  different  cases.  The  results  of  the  calculations  are  given  in  Appendix  B  and 
summarized  in  Table  1. 

It  is  interesting  to  note  at  this  point  that  the  techniques  of  orbit  reconstruction  discussed  in 
Sec.  Ill  satisfy  the  observability  test.  The  flight  times  between  points  in  orbit  depend  both  on 
radius  and  angle  information,  as  is  readily  verified.  Four  measurements  of  times  along  the  or¬ 
bit  are  sufficient  for  complete  observability. 

Some  recent  work2^  has  addressed  itself  to  the  problem  of  finding  optimal  solutions  (in  a 
least- square  sense)  to  ncar-circular  orbit  control.  For  optimum  control  theory  to  be  really 
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TABLE  I 

SUMMARY  OF  CONTROLLABILITY  AND  OBSERVABILITY  PROPERTIES 
OF  THE  LINEARIZED  ORBITAL  EQUATIONS 


Meosurements  per 

Orbit 

Continuous 

One 

Two 

Four 

Meosurement 

Observing  (in-plone) 

Rodiol  observotions 

No 

No 

No 

No 

Tongentiol  observotions 

No 

No 

Yes 

Yes 

Rodiol  ond  tongentiol 

No 

Yes 

Yes 

Yes 

Control  (in-plone) 

Rodiol  thrust 

No 

No 

No 

No 

Tongentiol  thrust 

No 

No 

Yes 

Yes 

Rodiol  ond  tongentiol 

No 

Yes 

Yes 

Yes 

Observing  (out-of-pl one) 

Orbit  normol  observotion 

No 

No 

Yes 

Yes 

Control  (out-of-plone) 

Orbit  normol  thrust 

No 

No 

Yes 

Yes 

useful,  the  system  must  be  completely  observable  and  controllable.  Thus  I.  B.  Chammas'  study 
deals  with  the  cases  which  are  completely  controllable  and  observable  (see  Table  1).  The  cost 
function  which  is  minimized  is 

('<* 

q  -  \  (u  ru  +  ay  'y)  dt 
vo 

where  y  is  a  linear  transformation  of  the  measured  vector  x*  The  author  studies  the  cases  in 
which  (1)  a  is  finite  and  large  and  a  =  0,  for  which  he  finds  a  time-varying  optimal  control,  and 
(2)  a  —  00  and  a  ^  0,  for  which  he  finds  a  constant  control  law. 

The  optimal  control  laws  depend  on  all  components  of  the  measurement  vector  x.  While  the 
instrumentation  presently  available  does  not  allow  us  to  conveniently  measure  nor  x^,,  it  is 
nonetheless  instructive  to  consider  the  optimal  control  policies  as  performance  bench-marks. 
Future  study  should  also  focus  on  performance  criteria  of  the  type 

q  -  1  (  |u|  +  ay'y)  dt 

Jo 

and  constraints  of  the  type 

pa 

\  |u|  dt  <  L 

%Jo 

C.  Stable  Control  of  the  System  with  Daily  Sampling 

It  is  important  now  to  examine  the  exact  meaning  of  the  lack  of  controllability  of  the  dynamic 
equations  in  the  case  of  one  sampling  per  orbit.  It  is  seen  from  ( IV- 23)  that  the  longitudinal 
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coordinate  is  controllable,  while  is  not.  The  question  therefore  arises  as  to  whether  it  is 
possible  for  x^  to  grow  without  bound  when  x^  only  is  controlled.  Physically,  this  would  corre¬ 
spond  to  the  orbit  becoming  increasingly  elliptical,  while  the  orbit  semimajor  axis  remains 
constant,  as  required  by  the  fact  that  x^  is  controlled.  We  shall  prove  that  the  double  thrusting 
policy  which  was  introduced  in  Sec.  Ill  has  the  property  that  the  orbital  eccentricity  remains 
constant  and  equal  to  its  initial  value. 

From  the  z-transfer  matrix  (IV-Z3),  the  transfer  function  between  x^  (the  longitudinal  co¬ 
ordinate)  and  u^  (tangential  thrusting)  is 


67TZ 


a;  (z  —  1 ) 


2  u2 


(IV-26) 


which  corresponds  to  a  double  integrator.  Thus  x^  is  controllable.  But  controllability  of  the 
longitudinal  coordinate  sampled  once  per  day  means  that  x^  can  be  driven  to  zero.  This  in  turn 
implies  that  the  orbital  period  and  therefore,  by  Kepler's  third  law,  the  semimajor  axis  of  the 
orbit  can  be  made  constant. 

It  is  easy  to  show  that  for  linear  expansions  of  the  orbital  motion  the  semimajor  axis  equals 
the  average  of  the  radius  over  one  orbit. t  Then  the  semimajor  axis  must  be  x^  +  r^  and  an  ex¬ 
pression  for  the  average  over  one  orbit  of  x  ^  is  found  from  the  transition  equation  [see  (IV- 10)], 

x .  =  (4  —  3  coswt)  x,(0)  +  i  sin  cot  x_  (0)  +  —  (1  —  coscot)x.(O)  ,  (IV- 27) 

1  1  CO  2  Cl?  4 


and  thus 


Xj  =  4x^0)  +  ^  x4(0) 

If  we  drive  x^  to  zero,  x^  must  equal  zero,  so  that  the  semimajor  axis  equals  r  . 
cessive  samplings  the  uncontrolled  coordinate  x^  is  constrained  by  the  equation 


(IV-  28 ) 
Then  on  suc- 


(I\  -29) 


Thus,  as  was  expected  on  physical  grounds,  the  values  of  x^  and  x^  are  not  independent  when  x^ 
is  controlled  to  zero. 

Consider  now  what  happens  if  one  samples  x^  so  that  its  value  equals  x^ .  Then  Eq.  (1V-28) 
becomes 


3x.  f  —  x„  =  0 
1  co  4 


(IV- 30) 


which  is  the  relationship  that  was  established  in  See.  II  when  sampling  was  performed  in  such  a 
manner  as  to  eliminate  the  daily  variation  in  r^  and  r^  [see  Eq.  (11-19)]. 

In  Sec.  Ill  we  had  established  that  the  daily  motion  (eccentricity  dependent)  was  eliminated 
if  one  observed  once  per  orbit  by  taking  the  average  of  two  observations  at  intervals  7r/to.  Also, 
it  was  claimed  that  splitting  the  thrust  value  (determined  once  per  orbit)  in  two  equal  parts  ap¬ 
plied  at  7T / Cl?  intervals  maintained  the  orbital  eccentricity  at  a  constant  value.  To  check  this,  we 
use  again  the  transition  matrix  (IV-  10)  and  find  immediately  that  the  new  observable 

*3  =  j  [*3(0)  +  x3(  ^  )]  (IV-  31  ) 


t  In  general,  the  semima|ar  axis  differs  fram  the  time  average  of  the  radius  taken  during  ane  orbit  (see,  for 
example,  Ref. 26). 
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does  not  contain  any  orbit  periodic  components.  This  new  observable  is  thus  the  mean  longitudinal 
position. 

To  find  the  effect  of  symmetric  thrusting  at  7r/u>  intervals  on  the  eccentricity  e,  we  note  that 


X 

_  max  o  1  max 
€  -  -  =  - 


1  min 


(IV-  32) 


Thus  for  the  eccentricity  to  be  constant  it  must  be  x.  =  constant. 

J  1  max 

From  the  transition  matrix  it  is  easily 


x,  =  4x .  ( 0 )  +  —  x.(0)  -f 

1  max  1  co  4 


x  (°)i2  - 

1  +  [5x^0)  t  ^  x4(0)J^ 


(IV-  3  3) 


Since  we  control  x^,  it  is  because  of  (1V-29), 


1  max 


rx  (0)iZ 

[4r  )  M3Xi(0)+  -  X^OH 


From  the  sampled  equations  for  T  =  7r/a>  (IV- 20)  used  with  the  constraint  u^fkT 
k  even,  one  has  easily 


T) 


(IV-  34  ) 
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(IV—  35) 


which  yields 


x2(kT  f  2T )  =  x2(kT) 


(IV-  36) 


and 


x .  ( k T  i  2T )  -  x  (kT)  +  -  u0(kT) 
1  1  u)  2 


x^  (kT  *  2T)  x4(kT)-6  u2(kT) 


(IV-  37) 
(IV- 38) 


Since  remains  unchanged  and  so  does  the  quantity 


3Xj  f  —  x^  (from  the  last  two  equations), 


the  eccentricity  remains  unchanged  because  of  (IV- 32)  and  (IV- 34).  Thus  we  have  shown  that, 
when  x^  is  controlled  by  thrusting  equally  at  points  180°  apart  in  orbit,  the  eccentricity  is  not 
changed.  Using  the  eccentricity  as  a  parameter,  if  the  sampling  is  done  so  that  0.  one  has 
from  (IV-  34),  (IV- 32),  and  (IV-  29), 

— e  v 


o 


x  A  -  2  co  e  r 
4  o 


(IV-  39) 

( I V- 40 ) 


for  this  control  policy. 
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To  conclude,  we  have  shown  how  it  is  possible  to  control  the  longitudinal  motion  of  the  satel¬ 
lite  by  sampling  once  per  orbit.  The  radius  is  uncontrollable,  but  constrained  to  the  value  of  the 
longitudinal  drift  rate  by  Eq.  (1Y-29).  Equation  (II- 9),  established  for  an  appropriate  value  of  the 
sampling  phase,  was  found  again  to  be  valid.  It  was  shown  how  it  is  possible  to  ensure  that  the 
eccentricity  remain  constant  and  thus  that  x^  (t)  and  x^t)  be  bounded,  =  -  e  r  ,  and  x^  2  to  c  r  . 

V.  ANALYSIS  AND  DESIGN  OF  A  FINITE-STATE  CONTROLLER 
A.  Preliminaries 

Section  V  describes  the  design  of  the  controller  which  was  actually  flown  in  LES-6.  The  de¬ 
sign  is  based  on  the  assumption  that  only  satellite  angular  position  can  be  measured  with  the  in¬ 
strumentation  available.  At  the  time  of  this  design,  it  was  assumed  that  daily  angular  drifts  of 
fractions  of  a  degree  per  day  could  not  be  measured  in  practice.  A  preliminary  estimate  of 
sensor  errors,  solar  clock  errors,  and  the  coupling  of  orientation  errors  into  position  errors 


led  to  t lie  conclusion  that  satellite  position  would 
2  degrees  peak- to- peak  over  long  periods  of  time 
established  that  longitude  would  be  measured  in 
window  was  set  at  2  degrees. 

It  was  also  decided,  on  the  basis  of  the 
analysis  of  Secs .  Ill  and  IV,  that  we  would  sam¬ 
ple  once  per  day  and  observe  only  the  mean 
daily  longitudinal  position.  Further,  it  was 
decided  that  the  daily  observation  would  con¬ 
sist  of  the  average  of  two  observations  made 
at  opposite  points  in  orbit  (with  respect  to  the 
earth  center  ),  to  decouple  the  observation  from 
the  daily  (eccentricity  dependent)  motion.  Fi¬ 
nally,  it  was  decided  to  thrust  symmetrically 
with  respect  to  the  earth  center,  so  that  ec¬ 
centricity,  although  not  controlled,  would  not 
increase,  as  discussed  in  Sec.  IV  (see  Fig.  5). 

Since  drift  velocity  was  not  directly  meas¬ 
urable.  damping  of  the  motion  could  be  a- 
chicved  by  observing  the  transitions  of  the  sal¬ 
t'll  ite  into  adjacent  two-degree- wide  longitude 
bands.  This  suggested  a  finite- state  machine 
realization  of  the  controller. 

The  thrusters  available  on  t ho  satellite 

ar  c'  of  two  types:  cold-gas  thrusters  and  pulsed 

plasma  thrusters.  The  gas  thrusters  provide  a 

-7 

and  the  plasma  thrusters  about  0.95  x  10  m/se 
linear-  acceleration  caused  by  the  earth's  equator 
Thus  it  was  decided  that  the  plasma  thrusters  wc 
their  maximum  effect  (in  continuous  operation)  i 


lot  be  measured  more  accurately  than  about 
with  an  rms  value  of  1 . 3  degrees.  It  was  thus 
quantized  manner  and  the  smallest  quantization 

SUN 

♦ 


PROJECTION  I 
OF  SUN  LINE  I 


Fig.  5.  Meosurement  of  sotellite  position  on  LES-6. 

-  5  2 

inear  acceleration  of  about  0.35  x  1()  m/sec 

2  i 

These  figures  arc  to  be  compared  to  the 

-  7  2 

ial  elliptic itv  which  is  at  most  0.545  x  io  m/sec 
uld  be  used  only  for  steady- state  tracking,  since 
only  about  twice  the  disturbing  acceleration  at 


t  Both  figures  keep  into  occount  the  duty-cycle  of  operarion  of  the  system. 


the  worst  earth  longitude.  The  need  to  operate  with  two  widely  different  thrusters  posed  an  ad¬ 
ditional  constraint  on  the  design  of  the  controller.  The  gas  thruster  system  was,  however,  con¬ 
sidered  the  primary  system,  since  there  has  been  no  previous  flight  experience  with  the  plasma 
thrusters.  As  discussed  in  Sec.  VIII,  these  were  indeed  extremely  successful  and  have  been  used 
exclusively,  up  to  the  time  of  this  writing,  for  steady- state  station  tracking. 

The  orbit  control  system  of  LES-6  has  the  following  modes  of  operation* 

Command  Mode:—  Either  the  gas  or  plasma  thrusters  can  be  operated  on  command 
from  the  ground.  The  earth  direction  and  the  firing  angles  (with  respect  to  the  orbit  tangent)  are 
determined  by  using  on-board  sensors.  The  thrusters  shut  off  automatically  after  predetermined 
time  intervals,  which  can  be  selected  from  the  ground. 

Automatic  Mode:—  Operation  is  completely  self-contained  with  either  thruster  system. 
Two  modes  are  available,  with  automatic  switchover:  (1)  new  station  acquisition,  and  (2)  station 
tracking,  with  self-contained  damping  of  the  relative  motion. 

In  all  modes,  thrust  is  applied  to  the  satellite  in  the  direction  of  the  orbit  tangent,  which  for 
a  circular  orbit  is  normal  to  the  satellite-earth  line,  determined  by  an  earth  sensor.  Total  im¬ 
pulses  of  different  values  are  obtained  by  firing  the  thrusters  for  a  fixed  part  of  each  of  several 
rotations,  the  spin  rate  being  measured  by  sun  sensors. 

The  satellite  position  is  determined  using  the  combination  of  sun  and  earth  sensors  described 
in  Sec.  111.  Sun  sensors  were  preferred  to  star  sensors  because  of  their  simplicity  and  rugged¬ 
ness.  The  position  of  the  satellite  must  then  be  compared  with  the  position  of  the  ground  station 
at  a  given  time.  Since  any  earth  meridian  moves  at  a  constant  angular  speed  in  space  (disregard¬ 
ing  small  secular  changes  of  earth  rotation  rate),  an  on-board  clock  can  determine  the  position 
of  a  station  of  known  longitude.  The  clock  delivers  a  pulse  every  12  hours:  the  time  difference 
between  a  clock  pulse  and  the  sensor  coincidence  pulse  is  then  a  direct  measurement  of  the  longi¬ 
tude  displacement.  The  fact  that  the  satellite  position  is  measured  with  respect  to  the  sun  pro¬ 
jection  on  the  orbit  plane  requires  that  the  clock  be  a  "true  solar  clock,"  which  should  indicate 

the  transit  time  of  the  sun  at  the  desired  earth 
longitude  (see  Appendix  A).  This  transit  time 
(ephemeris)  displacement  is  shown  in  l  ig.  6 
for  the  year  1969;  secular  changes  can  be  dis¬ 
regarded  for  satellites  having  a  lifetinle  of  a 
few  years.  If  the  spin  axis  is  normal  to  the 
orbit  and  the  orbit  is  equatorial,  the  true  solar 
time  of  the  earth  station  can  be  compaml  di¬ 
rectly  with  the  position  or  coincidence  pulse  de¬ 
livered  by  the  sensor  system.  A  high-precision 
oscillator  and  counter  deliver  a  pulse  every  mean  solar  day,  which  is  advanced  or  retarded  elec¬ 
tronically  in  accordance  with  the  day  of  the  yeart  (see  Sec.  Vll).  Alternatively,  a  "mean  solar 
clock"  can  be  used  if  the  position  of  the  satellite  is  measured  with  a  special  sun  sensor?  or  a  star 
tracker.  Spin  axis  misalignment  causes  a  position  error,  whose  worst-case  magnitude  during  a 


Fig.  6.  Sun  transit  time  displacement  for  1969. 


t  For  o  sotellite  lifetime  of  5  yeors,  clock  occurocy  of  obout  1  port  in  10  is  necessory. 
t  U.S.  Potent  Application  643493,  filed  5  June  1967. 
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year  is  ±0.55°  of  longitude  per  degree  of  misalignment  for  a  satellite  such  as  LES-6  (see  Ap¬ 
pendix  A  and  Kef.  31 ). 

Automatic  acquisition  of  new  stations,  i.e.,  new  earth  longitudes,  is  based  on  the  fact  that 
the  time  of  arrival  of  the  daily  pulse  from  the  internal  clock  determines  the  stationkeeping  longi¬ 
tude.  If  satellite  position  is  measured  at  the  9  a.m.  position  in  orbit  (point  B  in  Fig.  5),  the  con¬ 
trol  system  maintains  synchronism  with  that  earth  station  whose  true  solar  time  is  9  a.m.  at 
the  time  of  arrival  of  the  internal  clock  pulse.  Therefore,  new  stations  are  selected  automat¬ 
ically  by  changing  the  phase  of  the  internal  pulse,  i.e.,  by  resetting  the  clock  with  a  command 
when  the  true  solar  time  of  the  desired  station  is  9  a.m.  Station  changes  could  be  completely 
automatic  if  the  time  of  occurrence  of  the  reset  pulse  were  preprogrammed  on  the  satellite. 

After  acquisition,  the  satellite  can  stationkeep  at  any  geocentric  longitude,  which  implies 
bidirectional  control  action,  since  the  disturbing  force  changes  sign  in  accordance  with  Eq.  (11-36), 
and  its  motion  is  constrained  to  a  band  around  the  station,  whose  width  is  determined  by  the 
measurement  accuracy.  A  sequential  damping  procedure,  designed  to  use  fuel  with  high  effi¬ 
ciency,  ensures  the  convergence  to  this  band. 

The  most  important  performance  parameters  are  fuel  consumption  (both  in  steady  state  and 
in  damping),  the  settling  time  after  acquisition,  sensitivity  to  measurement  errors,  the  time 
necessary  to  acquire  a  new  station  and  the  number  of  possible  station  changes. 

B.  System  Analysis  (Ref.  32) 

In  this  section  we  use  the  "weak-coupiing"  approximation  that  was  established  in  Secs.  HI 
and  IV.  The  relevant  equations  are  repeated  here,  with  the  definition  of  constants  necessary  for 
the  applications.  Variables  are  sampled  once  per  day,  but  the  star  notation  is  dropped  for 
simplicity. 

The  geopotential  perturbation  is  given  by  ( i  1—  36), 

Ar ,  -  —  1.54X10  ^  sin  2  (A  A  )  deg/day^  ,  ( \  - 1  ) 

t  r  O 

where  A ^  -  73°  east  longitude  is  the  position  of  the  minor  axis  of  the  geopotential  ellipse.  The 
change  in  radius  is  related  to  the  angular  drift  velocity  by  (11-19), 

^6  2  Mtleg/day)  . 

r  3  360°  '  1  ' 

o 

where  r  synchronous  radius  ~  42, 165  km. 
o 

The  effect  of  the  thrusters  is  obtained  by  computing  the?  total  change  in  linear  velocity  per 
day  Av^,  knowing  the  thrust  level,  the  mass  of  the  satellite,  and  the  time  during  which  the 
thrusters  operate.  For  small  Av^,  the  daily  angular  acceleration  caused  by  the  thrusters  is  ob¬ 
tained  from  the  laws  of  orbital  motion, 

Av  .  7 

Arf  -  -3  X  360°  (deg/day  )  (V-3) 

o 

where  vq  =  synchronous  linear  velocity  ~  3,  300  m/sec.  The  total  angular  acceleration  is 
therefore 

A  A,,  +  A™  .  (V-4) 

u  1 
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A  convenient  approximate  representation  of  the  effect  of  the  geopotential  is  obtained  using 
a  "potential  energy  function"  V  ,  dependent  on  the  satellite  longitude  sampled  once  per  day.  By 
integration, 

V„(A)  -  \  X„dA  =  -4.  4  X  10"2  eos2(A  -  X  )  (deg/day)2  ,  (V-5) 

C  Jo  C  o 

recalling  the  necessary  transformation  from  radians  to  degrees.  The  thruster  is  designed  to 
produce  a  very  small  change  of  orbital  angular  velocity  per  day.  If  the  satellite  drift  per  day 
remains  small,  firing  the  thruster  every  day  for  a  fixed  time  results  in  a  constant  daily  change 
of  mean  drift  rate,  i.e.,  a  constant  mean  drift  acceleration.  If,  furthermore,  the  firing  time 
per  orbit  is  made  to  depend  only  on  the  satellite  mean  daily  longitude,  the  thruster- imparted  ac¬ 
celeration  is  sensibly  conservative,  and  a  second  "potential  energy  function"  can  describe 
the  longitude  caused  by  the  thruster.  The  satellite  drift  motion  can  therefore  be  studied  in  an 
approximate  but  expressive  way  with  the  aid  of  a  "total  potential  energy  function"  in  the  mean 
relative  satellite  longitude  X  measured  onee  per  day. 

In  the  system  discussed  here,  X  is  quantized  in  steps  of  magnitude  AL.  The  daily  thruster 
impulse  takes  on  two  different  values,  1*  between  AL  and  2AL  or  —  AL  and  —  2AL,  and  IM  else¬ 
where  for  X  |  <  L.  The  direction  of  thrust  is  such  as  to  eanse  the  satellite  to  drift  toward  the 
station.^  The  resulting  potential  energy  function  is  shown  qualitatively  in  Fig.  7. 


Fig.  7.  Quolitotive  behovior  of  the  potentiol  Fig.  8.  Switched  domping  technique, 

energy  function. 

If,  for  example,  the  satellite  were  initially  at  L  with  no  drift  velocity,  it  would  oscillate 
between  L  and  .  Ramping  is  therefore  necessary  and  is  accomplished  in  "damping  bands" 

A  and  A'  ,  AL  wide  about  the  station,  as  seen  in  Fig.  8.  When  the  satellite  is  in  A,  a  daily  im¬ 
pulse  is  applied  if  the  satellite  has  entered  A  by  crossing  the  zero  (station)  line,  and  no  impulse 
is  applied  if  it  has  entered  from  B.  A  similar  firing  strategy  is  used  when  the  satellite  is  in  A1. 


t  Notice  thot  on  increose  in  instontoneous  orbitol  velocity  increoses  the  semimojor  oxis  ond  the  period,  thus 
octuolly  slowing  down  the  sotellite  with  respect  to  the  stotion. 
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ACQUISITION  MODE 


Consider,  as  an,  example,  the  behavior  of  a  satellite  starting  at  P,  with.no  drift  velocity,  as- 

.1  '■  •  1  til  ■  ■  ')'•  ■  *  ■  !:  1.  up-ii  urfl  1o  nouoniR  t\  •  ‘7C-nq 

suming  for  simplicity  that  =  0.  The  satellite,  would  b?e  accelerated  by,Lhe  thyustcrs  and  reach 


the  station  with  kinetic  epergy  (in  the  variable^  A )  equaj  to  V^.  No  impuls^  \ypuld  ,be  applied  in 
A’,  but  upon  entering  A,  et  d*i}Tipi.ng  impulse  \yould  be  called  for  once  per  da/y,.  so  that ,  in.  bapd  A 


the  kinetic  energy  of  the  satellite  would  decrease  by  AV  (Fig.  8 /  ..^hus  the,  satellite  \yould  only 


reach  instead  of  invert  its  motion,  and  enter  A  from  B.  No  impulse  would  therefore  be 
applied  until  the  satellite  crossed  from  A  to  A\.  '  In  each^complete  oscillation/  the  energy  de¬ 
creases  by  2AV  until  the  satellite  remains  in  bands  A  or  A'.  As  seen  in  the  following,  the 
damping  impulse  is  changed  adaptively  to  reduce  the  residual  relative  kinetic  energy  to  very 
small  values.  I*  and  1M  are  made  large  enough  so  that  there  is  a  potential  energy  minimum  in 
a  neighborhood  of  the  station  for  any  station  longitude  and  a  stable  equilibrium  point  occurs  within 
±AL  from  the  station.  '  > 


1  Station  Acquisition  Mode 

The  automatic  system  enters  the  station  acquisition  mode  when  its  longitude  differs  from 
the  station's  by  more  than  ±L  degrees.  L  is  made  sufficiently  large  so  that  this  event  only  oc¬ 
curs  as  a  consequence  of  a  station  change.  Redundancy  techniques  are  used  to  avoid  spurious 
transitions  into  the  acquisition  mode.  The  thruster  immediately  fires  for  a  fixed  length  of  time 
at  opposite  points  of  the  orbit,  causing  the  satellite  to  drift  toward  the  station.  When  the  de¬ 
sired  geocentric  longitude  is  reached,  the  satellite  returns  to  near  synchronous  orbit  by  firing 
the  thruster  in  the  same  manner  in  the  opposite  direction.  Having  completed  the  homing  firing 
sequence,  the  system  switches  to  the  station-tracking  mode.  As  a  result  of  orbit  perturbations 
during  the  transfer  between  stations,  in  general  at  the  beginning  of  tracking  mode  operation  the 
satellite  has  a  residual  drift  velocity  relative  to  the  new  station,  which  must  be  damped. 

If  Aq  is  the  drift  rate  impressed  on  the 
satellite  at  the  beginning  of  the  acquisition 
mode,  then  within  the  limits  of  approximation 
of  ( V—  \  )  the  energy  conservation  theorems  in 
the  variables  A  and  A  can  be  written  as  for  a 
point  mass  whose  motion  were  described  by 
this  equation.  Prom  (V-5)  and  Pig.  9,  for 
acquisition  of  any  station  it  must  be  \  A^  > 

AV  -  -  8.8  x  10'2  (deg/day)2,  or  X  0.42 

deg/day.  Then  let  Aq  -  0.42  k  deg/day,  where 
k  >  1  is  a  coefficient  which  determines  the  ac- 
quisition  speed.  If  the  relative  drift  velocity 
imparted  by  the  thrusters  is  i  deg/day  per 
sec  of  thruster  operation,  the  thruster  on- 
time  for  acquisition  is  T^  =  0.42  k/i  sec.  The 
drift  velocity  must  vary  by  the  same  amount 
twice  because  there  must  be  two  equal  and  op¬ 
posite  impulses.  If  the  fuel  consumption  is 
m,  gm/sec,  the  fuel  required  per  station  acquisition  is  mf  =  0.84  m  k/i.  Notice  that  the  in- 

2  A 
equality  A  /2  : *  AV,*  need  be  satisfied  only  for  transfers  between  the  minor  and  the  ma  jor 

1  o'  ^  Ci  max  •  J 

axis  of  the  earth’s  equator.  One  could  compute  on  board  the  minimum  value  of  the  acquisition 


t  t 


MINOR  MAJOR 

AXIS  AXIS 

V _  J 

OF  EARTH'S  EQUATOR 

Fig.  9.  Kinetic  energy  for  the  transfer  between  orbits. 
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impulse  as  a  function  of  the  required  longitude  displacement,  using  (V-5),  and  fire  the  thruster 
accordingly,  if  the  additional  complexity  were  justified. 

The  speed  of  station  acquisition  ean  be  derived  by  integrating  (V-l )  with  initial  condition 
A  =  A^  for  various  longitudes  of  the  old  and  the  new  station.  To  obtain  results  in  closed  form, 
(V-5)  can  be  approximated,  dropping  the  phase  term  0qJ  by 

V(.(A)  s -4.4  x  10'2  ^1  -  t|-  X2j  (deg/day)2  (V-6) 

for  —  tf/4^  A  ^  tt/4,  where  A  must  be  expressed  in  radians.  This  second-order  approximation 
to  cos  2A  is  in  error  by  5  percent  at  A  =  it/8.  The  worst  ease  for  a  90°  transfer  is  clearly  that 
in  which  the  old  station  is  on  the  minor  axis  of  the  earth's  equatorial  ellipse  and  the  new  station 
is  on  the  major  axis  (Fig.  9).  The  initial  kinetic  energy  becomes 

V  =  8.8  X  1 0-2  k2  (deg/day)2  .  (V-7) 

2  2  2 

Defining  £  =  2k  —  1 6 ( A  /t t  ),  the  kinetie  energy  during  the  first  half  of  the  trajectory  is,  refer¬ 
ring  to  Fig. 9, 

V  (O)  =  4.4  X  1C )"2|  ,  (V-8) 

and  therefore  expressing  the  velocity  in  rad/day, 

i 

A  -  (0.297 7r/l 80)  ,  (V-9) 

The  time  to  travel  an  angle  of  45°  is  then 

,  or,  C*/4  _i 

tl  =  0.2977T  J  ^  2  dA  ( \  - 1 0 ) 

o 

which  integrates  to 


t,  -  152  sin'1  (Z~  z  k_1)  days  .  ( V-  4 1 ) 

Similarly,  in  the  second  half  of  the  trajectory, 

V'e(A)  =  4.4  X  10'2  [2k2  -  2  +  ^|  (X  -  |)2]  ,  ( V-  12  ) 

7 r 

and  defining 

>P  =  2Mk2-  l)l/[2k2-  1)1  -  1]  ,  ( V- 1 3 ) 

t^  152  log^  days  ,  (V-l 4) 

and  the  worst-case  time  for  a  90°  transfer  is 

tw(9()0)  =  tl  +  t2  days  .  (V-15) 


The  worst-ease  acquisition  time  for  transfers  smaller  than  90°  oeeurs  when  the  new  and  the  old 
station  are  loeated  symmetrically  with  respect  to  the  45°  longitude  (Fig.  9)  and  can  be  calculated 
easily  with  the  previous  procedure.  The  acquisition  time  for  180°  transfer  between  any  two  posi¬ 
tions  is  clearly  twice  tw  (90°).  The  best-ease  90°  transfer  occurs  when  the  old  station  is  on  the 
major  axis  and  the  new  station  is  on  the  minor  axis. 
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2.  Fuel  Consumption  in  Steady  State 

The  optimum  steady- state  condition,  from  the  point  of  view  of  fuel  consumption,  would  be 
one  in  which  a  constant  acceleration  were  applied  to  the  satellite  to  equal  the  tangential  drift 
acceleration  at  the  satellite  geocentric  longitude.  This  ideal  situation  cannot  be  achieved  because 
perfect  damping  is  not  possible  in  practical  controllers.  As  will  be  seen  later,  the  damping  im¬ 
pulse  is  reduced  adaptively  to  a  small  value  I  during  the  damping  maneuver.  After  the  controller 
has  reached  the  last  value  of  1  ,  either  the  impulse  in  A  and  A'  of  Fig.  8  is  reduced  to  zero  (first 
system),  or  the  impulse  remains  at  the  value  I  but  the  switching  is  discontinued  (second  system). 

a.  First  System 

No  impulse  is  applied  in  the  damping  band  after  damping  has  been  completed.  The  potential 
configuration  in  this  case  is  shown  in  Fig.  10.  The  longitude  interval  4AL  is  sufficiently  small 
so  that  the  external  acceleration  given  by  (V-l)  can  be  considered  constant.  Impulses  tl1  and  I* 
are  applied  in  opposite  directions  on  the  sides  of  the  damping  band.  For  a  daily  impulse  1  the 


Fig.  10.  Potential  energy  for  the  first  system. 


RELATIVE 

LONGITUDE 


2 

energy  acquired  by  the  satellite  in  traveling  AA  degrees  is  V  IAA  (deg/day)  .  Lot  F  be  the 

maximum  value  of  f  and  define  ot  _  1  /F  ,  and  a  '  I’/F  •  Since  the  satellite  must  remain 
e  n  n  e  e 

within  ±Ah  of  the  origin  after  damping,  a  1  must  be  greater  than  one.  At  the  end  of  the  damping, 
the  satellite  has  at  most  residual  energy 

VR=(In  +  fe)AL  '  (V'16) 

This  worst-ease  situation  oceurs  if  the  satellite  reaches  point  B  in  Fig.  10  with  zero  velocity. 

If  this  docs  not  oecur,  the  satellite  exits  from  the  damping  band  and  the  damping  cycle  continues. 
In  general,  the  satellite  velocity  will  be  zero  at  some  point  internal  to  the  damping  band,  in 
which  ease  the  residual  energy  is  smaller  than  V^.  The  damping  cycle  ends  at  the  origin,  since 
the  satellite  longitude  is  measured  in  steps  of  width  AL.  The  relative  kinetic  energy  is  zero  at 
p£,  where  the  potential  energy  is,  wfith  the  reference  of  Fig.  10, 
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Vc(P2)  =  VR  +  fcAL-(1n+2fe)AL  ' 


( V- 17 ) 


From  simple  geometric  considerations,  the  position  P 2  is  defined  by 

H"il>  I  .in  Ti]  i  »;!•:  ’if  •  J , ^  '  '  i  h  ■ '  i  '  '  t  •  1  j\  •  ■ 

x'  =  (1  +  2f  )  AI./U'  -  f  )  . 


n 

■  i  „■ 


e  1  e 

•i  Ut 


if  h  n* 


( V-  1 8 ) 


'Elementary  considerations  on  the  equation  of  motion  for  a  constant  acceleration  field  yield  that 
at  time 


ruc|  i  ■ 

■  / 


tA  /<v  xA/cr-  f  ) 


tv-19) 


where  A A  is  the  drift  velocity  with  which  the  satellite  passes  at  A  (Fig,  10), 
Pi. 

AA  =  [2(rn  +  2fe)ALF  . 


(V- 20) 


The  travel  time  between  A  and  B  is 
t/ 


lAB 


'*A  *B)//fe 


(V-21) 


where 


Ab  =  (21  AL?  ,  (V-22) 

The  velocity  at  is  zero,  and  therefore 

*B  V^V  •  (V-23) 

Sinee  this  eycle  repeats  itself,  the  fuel  consumption  per  second  is  proportional  to 

^  =  *lA  +  tB*^tA  +  lAB  +  lB  ^  ‘  (V-24) 

By  using  (V-17),  (V-18),  (V-22),  (V-24),  and  letting  A  -  XQ  -  6  in  (V-l),  the  expression  for  the 
steady- state  fuel  consumption  when  stationkeeping  at  any  position  in  orbit  is  obtained: 


Q 


(i) 


a'F 


Li  1  '2 


q  +  q + 


(cv  f  2  sin 26 
n  1  1 


(V-25) 


sin  26  | 


where  ~  (  nr  +  2  |  sin  26  |)  2  /(cv1  —  |  sin  26  |)  and 
for  sin  26  -*  0,  the  steady-state  fuel  consumption  for  f  0  is  found: 

,(i) 


Q  -  2 F  «  /(i  +  2  a  /a’) 
o  e  n7  n7 


/(a1  +  | sin  26  [ ) ,  By  taking  the  limit 


(V-26) 


b.  Second  System 

TIic  last  damping  impulse  1^  is  retained  in  the  part  of  the  damping  band  where  it  was  last 
applied.  The  potential  configuration  for  this  system  is  shown  in  Fig.  11.  The  residual  energy  is 
again  given  by  (V-l 6).  The  same  procedure  used  before  yields  the  steady- state  fuel  consumption 
for  this  system: 


q<2)  =  f 


g'q  *  %/q 


(a  +  2  I  sin  26  [  )2 

Z,  +  fo  +  — - j - ; - 

i  3  sin  26 


(V- 27 ) 
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Fig.  11.  Potential  energy  for  the  second  system. 


where  =  («n  +  |  sin  26  |  )2 ,  and  the  limit  for  sin26  -*  0  is 

Q  ^2  ^  -  2a:  F  /( 1 . 5  +  a  /a') 
m  n  c7  tv 


RELATIVE 

LONGITUDE 


(V-28) 


c.  Comparison  of  the  Two  Systems 

Comparing  Fqs.  ( V-2  6 )  and  (V-28 )  yields  Q^2  ^  ^  for  a  <  ^  cv  1 .  It  has  been  seen  that 

rv1  must  exceed  unity  to  ensure  the  existence  of  a  potential  minimum  in  a  neighborhood  of  the 

station.  Also,  must  be  less  than  unity  to  make  the  residual  energy  small.  In  general, 

this  inequality  is  satisfied  and  the  second  system  is  preferable  when  f  =  0.  For  instance,  for 

a'  =  2  and  a  -  1  ,  Q^'  =  0.4  F  and  Q;^2^  =  0.308  F  .  The  second  system  is  superior  when 
n  4  °  e  o  e  12)  (i) 

the  external  force  is  zero  or  sufficiently  small.  As  -*■  0,  /Q^  \  ■  When  the  external 

force  is  maximum,  i.e.,  I  sin  26  I  =  1,  then  =  1.27  F  and  =  1.16  F  .  Therefore,  the 

seeond  system  consumes  less  fuel  also  when  the  force  is  maximum.  Notice  that  when  0, 

and  tend  to  f  ,  i.  e.,  the  systems  tend  to  a  minimum  fuel  system  in  which  the  impulse 

per  orbit  applied  by  the  motor  equals  the  acceleration  per  orbit  due  to  the  external  force.  It  is 

therefore  possible  to  conclude  that: 

(1)  The  second  system  is  superior. 

(2)  The  impulse  V  =  a  1  Fe  should  be  the  smallest  compatible  with  the 
requirement  of  being  sufficiently  greater  than  Fe  to  give  some  safety 
margin.  A  good  value  is  approximately  c*  1  =  2,  the  final  selection 
being  made  on  the  basis  of  convenience  in  the  electronic  implementation. 

(3)  The  impulse  ln  in  the  last  step  of  the  damping  cycle  should  be  as  small 
as  possible. 

3.  Active  Damping 

After  the  satellite  has  acquired  a  station,  its  kinetic  energy  is  decreased  by  an  active 
damping  system.  A  first  requirement  is  to  have  a  high  potential  wall  (Fig.  7)  which  makes  it 
impossible  for  the  satellite  to  return  to  the  station  acquisition  mode.  This  determines  the  width 
L  of  the  tracking  band  and  the  impulse  level  between  +AL  and  —  F,  since  it  must  be,  from  (V-5), 
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^.0.088  (deg/day)  ,  where  Vj.  is  the  potential  energy  at  a  longitude  ±L  from  the  station. 
Letting  h  be  the  number  of  longitude  subdivisions  of  width  AL  in  which  the  applied  impulse  has 
the  value  I",  and  defining  a  M  similarly  to  a', 

V.  =  (a '  +  a  "  h)  F  AL  ^  0.088  (deg/day)2  .  (V-29) 

L  e 

Since  AL  is  chosen  on  the  basis  of  measurement  errors  and  of  system  requirements,  and  a1  was 
obtained  previously,  Eq.(V-29)  determines  the  product  a"  h. 

To  obtain  a  small  residual  value  of  the  drift  velocity  at  the  end  of  the  damping  process,  it 
is  necessary  either  to  measure  the  satellite  drift  velocity  accurately  or  to  reduce  the  damping 
impulse  adaptively.  In  the  first  case,  the  residual  drift  velocity  is  determined  by  the  measure¬ 
ment  accuracy;  in  the  second,  by  the  smallest  value  of  damping  force  which  can  be  used  in 
practiee.t  The  latter  approach,  which  leads  to  a  sequential  design  of  the  controller,  was  taken 
because  of  the  difficulty  of  measuring  drift  velocities. 

Two  damping  systems  will  be  compared,  in  which  successive  values  of  the  damping  impulse 
obey  a  relationship  of  the  type 

I.  =  ±1  a_k  ,  a  >  \  ,  (V-30) 

k  o 

and  the  sign  is  chosen  so  as  to  oppose  the  drift  velocity.  The  two  systems  differ  in  the  manner 

in  which  the  transition  from  1.  to  1,  ,  occurs: 

k  k  +  1 

(a)  Zero  Crossing  Adaptation:  The  transition  from  1^  to  lp.^  occurs  at  eaeh 
zero  crossing,  i.e.,  each  time  the  measured  position  of  the  satellite 
passes  from  region  A  to  A'  in  Fig.  8,  or  vice  versa. 

(b)  Convergence  Adaptation  The  transition  from  1^  to  j  occurs  eaeh  time 
the  satellite  converges  to  the  band  ±  AL  around  the  station.  In  general, 
several  zero  crossings  can  occur  before  convergence. 

In  studying  both  systems,  the  initial  kinetic  energy  will  be  taken  to  equal  Vrj  from  (V-29)  for 
a'  -  2,  a"  =  6,  and  h  =  6.  The  following  analysis  is  therefore  a  worst-ease  study  of  convergcnee 
time.  Two  definitions  are  useful  in  the  study  of  the  active  damping  system.  The  convergence 
time  is  the  time  after  which  the  total  energy  (in  the  X  variable)  Vj  remains  below  a  small  value 
whieh  is  the  residual  energy  toward  which  the  system  is  designed  to  converge.  The  damping 
efficiency  is  the  ratio  of  the  initial  kinetic  energy  (in  the  X  variable)  of  the  satellite  to  the  energy 
(in  the  X  variable)  expended  by  the  thruster  in  getting  the  total  energy  below  V^. 

Zero  Crossing  Adaptation:—  In  this  system,  the  damping  impulse  is  divided  by  a 
faetor  a  >  4  every  time  the  satellite  passes  through  zero.  After  k  half  oscillations  of  the  sys¬ 
tem,  the  total  energy  is 

k-  4 

vk  =  vL-«0FeAL  £  a_1  •  k^i  >  (V-31) 

i=0 

where  a?  F  is  the  initial  damping  impulse.  The  sequence  is  ended  when 

V  <  a  a-n  F  AL  ,  (V-  32) 

n  o  c  7 


t  Energy  is  coupled  into  the  system  by  errors  in  determining  the  boundories  of  longitude  bonds,  os  discussed  in 
the  following. 
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or 


n-  1 

V,  -  a  F  AL  T.  a" 1  <  a  a"nF  AL  .  ( V-  33 ) 

L  o  e  u  o  e 

i=0 

Using  ( V-29 )  and  executing  the  summation  in  (V-33),  rearranging  and  letting  a~n  =  €,  the  condi¬ 
tion  for  convergence  is 

aQ  >  38(a  -  1  )/(a  -  c)  .  (V-34) 

When  this  condition  is  not  satisfied,  the  rate  of  decrease  of  the  damping  impulse  is  too  great  in 
relation  to  that  at  whieh  energy  is  subtracted  from  the  system,  and  convergence  does  not  occur. 

If  ean  be  increased  sufficiently,  inequality  (V-34)  can  always  be  satisfied.  However,  the 
initial  damping  force  cannot  be  arbitrarily  large.  The  operation  of  the  adaptive  damping  sys¬ 
tem  will  eventually  be  limited  by  the  accuracy  of  the  longitude  measurement,  since  a  longitude 
error  6 L  in  establishing  the  boundary  of  the  damping  band  eauses  a  variation  of  the  potential 
energy  of  the  system 

6  V  -  o  F  SL  •  ( V-  35 ) 

o  e 

Requiring  that  6V  be  smaller  than  some  percentage  i}  of  Vj  ,  an  upper  bound  for  o>o  is  determined; 
e.g.,  for  5  L  =  0.25  °  and  tj  2%,  6.  Using  this  condition  in  ( V- 34 ),  a  ^  1. 19  —  6c/32.  Ob¬ 

viously,  operating  with  these  low  values  of  a  is  not  very  efficient.  The  energy  expended  by  the 
active  damper  on  the  satellite  is 

n 

W  =  2  V.  +  2V,  +...  +  2V  =  2V.  +2  Y.  V.  ( V-  36 ) 

L  j  l  n  i  j  k 

k  =  l 


or,  from  (V-  31 ), 


W  2(n  +  1)  V 


L 


2a 


a(n  4  1 ) 
a  —  1 


F  AL  +  2a 
e  o 


a  —  1 


)  (1 


a-(n+1))  F  AL 
c 


(V- 37) 


For  example,  if  e  l/32,  which  is  a  reasonable  value  to  have  a  sufficiently  small  I  ,  n  ~  19. 

With  the  values  n  =  19,  aQ  6,  =  38  F  AI,,  and  a  =  1.18,  the  damping  efficiency,  from  (V-37), 

is  r 7  =  0.08  5. 

Convergence  Adaptation:—  In  this  type  of  damping  controller,  is  retained  until 
the  satellite  is  trapped  in  the  damping  band  ±AL  wide  around  the  origin.  In  the  rest  of  this  sec¬ 
tion,  the  most  efficient  damping  sequence  is  determined  for  dampers  of  the  second  type.  Con¬ 
sider  a  damping  sequence  in  which  =  V^a  ^  a  ^  convergence  is  obtained  after 

Nq  half-eyelcs  of  satellite  position.  From  the  previous  definitions  and  for  -  38  F^  AL,  the 
energy  expended  during  this  first  phase  is 

W  2V,  +  2f V ,  -  V  )  +  2f V  -  2V  )+....  +2  [VT  -  (N  -  1)  V  ) 

L  L,o  L  o  L  o  o' 

-  2N  V.  —  N  (N  -  1)  V  .  ( V-  38 ) 

o  L  o  o  o 


After  the  first  convergence,  the  number  N  of  half-eyeles  of  satellite  position  necessary  to  obtain 
eacli  sueeessive  partial  convergence  obviously  is  the  smallest  integer  >  (a  —  1).  The  energy  W" 
expended  is  then,  for  a  >2, 
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Fig.  12.  Simplified  block  diogrom  of  the  Automotic 
Orbit  Control  system  for  LES-6. 
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Fig.  13.  Block  diogrom  of  longitude 
meosurement  logic. 


Fig.  14.  Flow  diogrom  of  longitude 
meosurement  logic. 
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(V-39) 


n-  1 

W"  =  2N  Tj  V.  =  2Na(a-  l)'1  (1 -€  )  VQ 
i  =  0 

The  lowest  value  of  W"  is  obtained  for  a  =  2.  Sinee  the  energy  W'  is  independent  of  a,  the 

optimum  is  obtained  when  a  =  2.  Using  (V-38),  (V-39),  and  recalling  that  Vj  =  38  F^AL,  VQ  = 

cv  F  AL,  and  N  =  38 /a  , 
o  e  o  o 


* L) 


+  38  +  2(1  -  €)  cvo 


Na 
a  —  1 


(V-40) 


whieh  yields  =0.13  for  =  6,  a  =  2,  e  =  l/32,  and  N  =  1.  Notice  from  (V-39)  that,  if  1Q 
could  be  made  larger,  so  that  VQ  =  V^,  the  damping  efficiency  in  the  optimum  configuration 
would  be 


1 


'D  max  2(1  —  € )  Na 


0.25 


(V-41) 


This  efficiency  is  also  obtained  whenever  the  initial  energy  is  small,  so  that  the  system  always 

-  k 

converges  in  one  cycle.  The  values  of  the  successive  damping  impulses  are  then  =  —  2  I  . 
Sinee  cvo  is  selected  on  the  basis  of  the  maximum  potential  variation  caused  by  measurement 
errors,  the  only  parameter  to  be  determined  is  the  value  e,  i.e.,  the  final  value  1  ,  which  should 
be  retained  after  convergence  to  ensure  minimum  steady- state  fuel  consumption  and  should  be 
as  small  as  is  compatible  with  practical  requirements. 


VI.  REALIZATION,  SIMULATION,  AND  TESTING  OF  FLIGHT  SYSTEM 

A.  Logic  Design  of  Controller 

The  over-all  bloek  diagram  of  the  stationkeeping  system  is  shown  in  Fig.  12.  The  solar 
cloek  determines  the  position  of  the  station  whieh  is  selected  by  command  from  ground.  Sun 
and  earth  sensors  and  the  associated  coincidence  eircuitry  determine  the  satellite  position.  The 
longitude  difference  in  sign  and  absolute  value  is  stored  in  a  memory.  Measurements  are  made 
twice  per  orbit  at  satellite  anomalies  180°  apart  and  averaged.  The  controller  remembers  the 
sequence  of  longitude  difference  measurements  and  determines  the  required  thrust  level  in  ae- 
eordanee  with  the  laws  established  in  See.  V.  The  thruster  firing  logie  turns  the  thruster  on  so 
that  the  resulting  thrust  vector  is  always  tangential  to  the  orbit  and  has  an  appropriate  direction. 

A  command  from  ground  can  turn  the  thruster  on  at  any  time  for  the  initial  orbit  adjustments 
(See.  111-E).  Commands  ean  also  be  sent  to  the  cloek  to  initiate  it  to  the  appropriate  clay  of  the 
year(see  Appendix  A).  The  system  weight  is  about  4  pounds  and  the  power  consumption  is  120  mW. 


1.  Longitude  Difference  Measurement 

The  bloek  diagram  of  the  longitude  difference  measurement  logic  is  shown  in  Fig.  13.  The 
inputs  are  the  station  and  the  satellite  position  signals  E^  and  E^  and  P^  and  P  respectively. 

The  purpose  of  the  state  memory  is  to  discriminate  against  wrong  sequences  such  as  E^,  P^, 
Ej.  These  sequences  can  oceur  at  the  beginning  of  the  stationkeeping  experiment  and  whenever 
a  station  is  changed.  The  system  will  wait  until  P^  or  E^  followed  by  or  P  is  received  be¬ 
fore  enabling  the  longitude  counter. 

The  state  memory  also  determines  the  sign  of  the  longitude  difference  measurement  by  adding 
algebraically  the  difference  measurements  P^  and  P^,  as  results  from  the  flow  diagram  of  Fig.  14. 
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In  the  Early  and  Late  states  the  counter  counts  up.  In  the  Contrary  1  state  the  counter  counts 
down  until  it  crosses  zero,  then  reverses  direction  and  complements  the  sign  register.  The 
counter  has  an  overflow  register  which  determines  when  the  satellite  is  out  of  the  tracking  band. 

The  longitude  difference  measurement  is  transferred  in  fixed  increments  to  a  secondary 
memory  upon  receipt  of  a  transfer  pulse  which  is  generated  after  the  second  measurement  and 
the  thrusting  have  been  completed.  The  secondary  memory  contains  the  information  "the  satel¬ 
lite  is  {or  is  not)  in  band  A,  B,  C,  or  in  the  acquisition  band."  Signal  redundancy  techniques  are 
used  to  avoid  a  single  spurious  pulse  bringing  the  controller  into  the  Acquisition  Mode. 

2.  Control  Unit 

The  block  diagram  of  the  control  unit  is  shown  in  Fig.  15.  Upon  arrival  of  an  Acquisition 
signal,  a  pulse  train  is  fed  to  the  firing  timer,  which  opens  after  the  thrust  logic  has  sent  an 
Open  signal.  The  input  P£  is  used  to  count  days,  since  this  thrust  must  be  applied  for  two  con¬ 
secutive  days.  The  opposite  thrust  is  applied  when  the  satellite  arrives  at  the  damping  band. 

For  redundancy  and  to  avoid  false  signals,  the  logic  waits  for  two  B  signals,  or  a  transition 
from  Early  to  Late,  or  vice  versa.  When  either  situation  occurs,  the  opposite  thrust  is  applied 
for  two  days;  when  the  thrusting  operation  is  completed,  the  satellite  returns  to  the  station 
tracking  mode. 

In  the  tracking  mode,  the  operation  of  the  controller  is  determined  by  the  sequences  of  the 
inputs  A,  B,  C,  together  with  their  "sign"  E,  L.  The  "memory- state  diagram"  for  the  control¬ 
ler  in  the  tracking  mode  is  shown  in  Fig.  16.  The  thrust  output  is  determined  by  the  memory 
state  and  by  the  current  value  of  the  quantized  longitude;  thus  several  "output  states"  are  con¬ 
densed  in  one  "memory  state"  for  simplicity  of  presentation;  the  complete  output  state  diagram 
consists  of  6l  states.  From  the  acquisition  mode  the  controller  goes  to  the  state  (Damp  —  )  if 
it  is  early  and  to  (Damp  +  )  if  it  is  late.  In  these  states  no  damping  force  is  applied  whenever  the 
satellite  is  in  bands  +A  and— A,  respectively.  Whenever  there  is  a  partial  convergence,  the 


Fig.  15.  Block  diogrom  of  the  controller.  Fig.  16.  Memory-stote  diogram  of  the  controller. 
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sequenee  —  A,  +A,  —  A  occurs  if  the  controller  is  in  (Damp—)  and  the  sequence  +A,  —  A,  +A  if  it  is 
in  (Damp  +  ).  On  the  transition  from  the  states  Converge  to  the  states  Wait,  the  damping  foree 
is  decreased  unless  it  has  already  reached  its  minimum  value.  The  controller  exits  from  the 
states  Converge  upon  receipt  of  a  B  signal,  and  the  damping  force  returns  to  its  maximum  value 
upon  receipt  of  the  signal  Begin,  whieh  happens  after  two  successive  C  signals. 

The  adaptation  memory  input  is  the  Decrease  signal  (Fig.  15).  The  damping  foree  is  obtained 
by  ehoosing  a  frequency  in  the  frequency  selection  matrix,  whieh  is  sent  to  the  fixed  length  firing 
timer.  When  the  satellite  is  in  bands  ±B  or  ±C,  the  foree  level  is  established  in  a  similar  man¬ 
ner.  An  Inhibit  signal  prevents  the  operation  of  the  tracking- eontrol  mode  when  the  satellite  is 
in  Acquisition. 

When  both  the  firing  sequenee  and  the  longitude  difference  measurement  have  been  completed, 
a  Reset  and  a  Transfer  pulse  are  generated  and  fed  to  the  auxiliary  longitude  difference  memory 
in  Fig.  1  3. 


start  thrust 

STOP  THRUST 


3.  Thrust  Logic 

The  thrust  logie  block  diagram  is  shown 
in  Fig.  17.  The  earth  center  detector  delivers 
a  pulse  whenever  a  fixed  marker  point  on  the 
satellite  periphery  is  eollinear  with  the  earth 
and  satellite  eenters.  Sinee  the  firing  angle  of 
the  thruster  is  18°,  the  thruster  line  of  action 
is  plaecd  81  °  away  from  the  marker  point. 

After  the  arrival  of  the  earth  eenter  pulse,  if 
the  satellite  is  early  the  thruster  is  turned  on 
and  remains  on  until  the  firing  angle  eounter 
overflows.  The  sun  clock  used  to  measure  the 
firing  angle  delivers  512  pulses  per  satellite 
spin.  If  the  satellite  is  late,  the  thrust  must 
oppose  the  satellite  velocity  veetor.  A  delay 
of  180°  of  satellite  rotation  is  therefore  intro¬ 
duced  before  the  thruster  is  turned  on. 

The  logie  for  the  Plasma  Thrusters  differs  from  this,  in  that  there  are  four  thrusters  lo- 
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cated  90°  away  from  eaeh  other.  The  sun  eloek  is  used  to  gate  the  appropriate  thruster,  and 
the  180°  delay  selects  the  direction  of  firing. 


Fig.  17.  Thrust  logic  black  diagram. 


4.  Electronic  Solar  Clock 

In  order  to  compare  the  satellite  and  station  positions  in  orbit,  it  is  necessary  to  have  a  solar 
eloek,  as  discussed  in  See.  Ill. 

The  difference  between  the  mean  and  true  solar  transit  times  is  shown  for  the  year  1969  by 
the  smooth  curve  in  Fig.  18.  It  may  be  assumed  that  this  eurve  is  valid  for  subsequent  years, 
sinee  the  seeular  variation  of  the  difference  is  small  during  the  satellite  lifetime.  The  clock  is 
a  sequential  network  whieh  delays  the  output  of  a  eounter  by  a  time  D  (Fig.  18)  as  shown  in  the 
bloek  diagram  of  Fig.  19-  Two  pulses  per  mean  day,  12  hours  apart,  start  a  delay  eounter  whieh 
counts  minutes  until  its  contents  equal  the  contents  of  the  accumulator.  At  this  time,  the  output 
pulses  or  S2  are  generated  and  the  delay  eounter  is  reset  and  opened  until  the  next  pulse  from 
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Fig.  18.  Approximation  to  the  sun  transit  time.  Fig.  19.  Block  diagram  of  the  electronic  solar  clock. 

the  main  counter  arrives.  The  accumulator  contains  the  amount  of  delay  measured  from  the 
reference  line  in  Fig.  18.  As  shown,  the  reference  line  differs  from  the  mean  solar  day  line  by 
17  minutes.  This  is  taken  into  account  by  sending  the  alarm  command  at  an  appropriate  time. 

The  value  stored  in  the  accumulator,  which  is  an  up-down  counter,  is  computed  by  the  sequential 
network  shown  in  Fig.  19.  The  staircase  curve  in  Fig.  18  shows  the  actual  correction  obtained 
with  this  logic.  The  approximation  has  been  selected  by  a  trial-and- error  procedure,  with  a 
view  to  simplifying  the  logic  design.  Since  the  information  contained  in  the  solar  clock  logic 
must  be  preserved  when  the  satellite  is  in  shadow,  and  the  main  counter  must  continue  to  oper¬ 
ate  throughout  the  shadow  period,  the  system  is  powered  from  the  solar  bus  or  a  battery. 

The  most  important  function  performed  by  command  from  ground  is  the  selection  of  the  sta¬ 
tion.  This  is  done  by  resetting  the  main  counter  at  9  a.m.  local  mean  solar  time  of  the  station 
to  be  tracked,  since  the  satellite  position  is  measured  at  9  a.m.  satellite  true  solar  time.  Since 
the  daily  correction  is  obtained  as  a  delay  from  a  reference  line  17  minutes  earlier  than  the 
mean  solar  day,  the  time  at  which  the  main  counter  must  be  reset  is  8:43  a.m.  station  mean 
solar  time.  The  longitude  of  the  station,  the  longitude  of  the  command  site,  and  the  propagation 
delays  are  known  and  therefore  determine  the  site  mean  solar  time  at  which  the  command,  which 
is  called  ALARM,  must  be  sent.  An  error  in  the  time  at  which  the  ALARM  is  sent  only  causes 
an  offset  error  in  the  operation  of  the  control  system.  The  magnitude  of  this  offset  is  obviously 
1°  of  longitude  per  4  minutes  of  time  error,  so  that  very  reasonable  accuracies  can  be  achieved. 
The  command  system  of  the  solar  clock  also  performs  the  auxiliary  function  of  choosing  the  initial 
conditions  for  the  sequential  network  and  for  the  accumulator.  This  is  necessary  because  the 
actual  date  of  launch  is  not  exactly  known  beforehand  and  is  also  useful  as  an  additional  element 
of  flexibility  in  the  system. 

B.  Computer  Simulation  and  Flight  System  Testing 

An  extensive  computer  simulation  study  was  carried  out  during  the  program,^  to  yield  es¬ 
timates  (1)  of  the  fuel  consumption  when  operating  with  the  cold  gas  thrusters  in  the  presence  of 
random  measurement  errors,  and  (2)  of  the  dynamic  behavior.  In  the  course  of  this  work,  it 
was  found  that  it  was  advantageous  to  increase  the  width  of  the  B  longitude  band  from  2°  to  4°, 
since  this  resulted  in  lower  fuel  consumption  in  the  presence  of  random  errors.  This  modifica¬ 
tion  was  therefore  incorporated  in  the  hardware  design. 

The  Stationkeeping  Simulation  program  performs  the  following  operations  for  each  orbit: 

(1)  Integration  of  the  equation  of  motion  with  earth  tesseral  harmonics  through 
third  order. 
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(2)  Generation  of  the  coincidence  pulses  at  the  appropriate  sun  ephemeris 
time, 

(3)  Determination  of  the  satellite  position,  using  the  reconstruction  of  the 
sun  ephemeris  implemented  by  the  on-board  logic,  or  by  the  special 
sensor, 

(4)  Adding  random  error  {uniformly  distributed,  uncorrelated,  ±0.5°)  to 
the  measurement, 

(5)  Determination  of  the  thrust  level  required  for  the  next  orbit,  using  a 
subroutine  which  duplicates  the  behavior  of  the  finite-state  controller. 

The  outputs  are  the  osculating  orbital  elements  for  each  orbit  and  the  fuel  consumption.  It 
is  possible  to  simulate  operation  with  both  the  Plasma  and  Cold  Gas  Thrusters.  This  simulation 
program  can  take  as  inputs  the  orbital  elements  determined  by  radar  tracking. 

A  variable  step-size  Hamming’s  integration  routine  was  used.  It  is  important  to  use  a  var¬ 
iable  step- size  method  because  of  the  discontinuous  nature  of  the  forces  involved.  Much  work 
was  done  to  study  the  effect  of  integration  errors.t  This  is  most  important,  since  the  number  of 
successive  integrations  is  quite  large  when  the  simulation  is  carried  on  for  the  equivalent  of  five 
years  of  orbital  life.  In  order  to  handle  discontinuous  thrust  properly,  one  would  like  to  use  a 
small  step  size;  this,  however,  increases  the  number  of  necessary  integration  steps.  An  ap¬ 
propriate  compromise  was  found  and  its  validity  checked  by  comparing  the  results  given  by 
Hamming's  and  Kunge-Kutta  variable  step-size  routines. 

An  extensive  analysis  of  fuel  consumption  was  made  for  different  values  of  the  stationkceping 
longitude.  Uniformly  distributed  uneorrelated  noise  was  used,  since  it  seemed  a  sufficiently 
severe  test  for  the  system. 


Fig.  20.  Simulation  results  at  the  worst  longitude.  Fig.  21.  Simulation  results  at  the  unstable  point. 

The  design  constraint  of  not  being  able  to  measure  the  angular  drift  except  in  very  largely 
quantized  longitude  bands  mades  the  system  rather  sensitive  to  fast,  uneorrelated  noise.  The 
reason  for  this  is  that  the  controller  can  go  very  rapidly  through  the  damping  eycle,  as  a  conse¬ 
quence  of  spurious  sequences  of  signals  like  +A,  — A,  +A  or  — A,  +A,  —  A.  As  is  seen  in  See.  V,  the 
controller  responds  to  this  sequence  by  switching  to  a  lower  level  damping  force.  Figures  2  0 
and  21  show  typical  computer  runs  for  stationkceping  at  the  worst  longitude  (-295°  East)  and  at 
the  unstable  point  (-259°  East)  with  the  Cold  Gas  Thrusters,  with  random  error.  In  most  cases, 
the  longitudinal  excursion  is  within  ±3  degrees  of  the  station. 

Several  runs  were  made  to  estimate  the  fuel  consumption  in  the  presence  of  random  errors. 
As  is  to  be  expected,  the  fuel  consumption  is  a  function  of  the  stationkeeping  longitude.  Typically, 


t  E .  H.  Swenson,  unpublished  report. 
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the  system  will  use  2.6  kg  when  stationkeeping  for  five  years  near  a  point  of  maximum  tangential 
force,  and  0.6kg  near  the  unstable  point. 

Testing  of  the  system  was  carried  out  by  exciting  the  earth  and  sun  sensors  using  a  special- 
purpose  Simulator.  This  device  simulated  sequences  of  sun  and  earth  pulses  at  the  spin  rate  of 
the  satellite.  Circular  orbital  motion  was  simulated  by  ehanging  the  sun- earth  phase  uniformly 
in  time.  The  Simulator  had  provisions  to  change  the  orbital  speed  and  to  "jump"  to  any  point  in 
orbit  using  pushbuttons.  The  stationkeeping  system  clock  could  be  sped  up,  so  that  the  Day 
Pulses  would  occur  with  variable  repetition  rates,  down  to  as  little  as  1  minute.  In  this  manner, 
by  measuring  the  time  elapsing  between  the  Day  Pulses  E^,  E^,  which  were  available  externally 
via  telemetry,  the  "satellite  position"  on  the  simulator  could  be  changed  so  that  the  sensor  co¬ 
incidences  P^,  would  occur  at  the  appropriate  times  with  respect  to  the  clock.  By  changing 
the  time  between  E^  and  P^  (E^,  P^ )  for  a  given  clock  speed,  the  system  measures  different  lon¬ 
gitudes.  This  type  of  checkout  was  conducted  also  during  thermal- vacuum  tests. 

A  seeond  technique,  which  was  most  useful  in  testing  the  integration  of  the  sensors  and  logic 
systems,  was  employed  in  which  the  satellite  was  actually  spinning.  The  light  from  a  lamp  sun 
simulator  (appropriately  folded  by  mirrors)  was  used  to  excite  the  sun  sensors.  The  earth  image 
was  projected  on  screens  in  sueh  a  manner  as  to  obtain  approximately  the  right  earth  albedo. 

The  earth  image  could  be  moved  (discontinuously)  so  that  sensor  coincidences  could  be  obtained 
at  the  right  time  with  respeet  to  the  clock  pulses  also  in  this  ease. 


VII.  EARLY  FLIGHT  RESULTS  AND  COINCIDENCE  ANALYSIS 

The  longitude  history  of  the  satellite  from  26  September  1968  until  25  dune  1969  is  shown 
in  Fig.  22.  The  Automatic  Stationkeeping  System  operated  from  22  December  1968  to  the  end  of 
the  period.  Telemetry  was  recorded  daily  from  the  Automatic  System  longitude  sensors.  The 
ephemeris  transit  time  has  been  subtracted  from  the  sensor  data. 


Fig.  22.  Longitude  history  of  LES-6  since 
launch,  as  measured  by  the  an-baard  equip¬ 
ment  daily  (smoothed  data). 


Orbit  determinations  were  made  (independently  of  the  sensor  data)  onee  a  month,  and  more 
frequently  for  calibration  purposes  at  the  beginning  of  the  experiment.  These  are  indicated  by 
the  open  circles  in  Fig.  22. 

During  the  period  from  launeh  to  22  December  1968,  the  system  operated  mostly  in  the 
"Ground  Command"  mode.  During  this  time  the  following  operations  were  performed:  (a)  initial 
trimming  of  the  orbit,  using  the  Cold  Gas  Thrusters,  (b)  experimentation  with  the  Plasma  Thrust¬ 
ers,  to  determine  their  thrust  level,  and  (c)  experimentation  with  the  Automatic  System,  closing 
the  loop  for  brief  periods  of  time  and  monitoring  its  behavior  continuously. 
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All  operations  were  performed  successfully.  The  initial  orbit  (semimajor  axis  42,068  km; 
eeeentrieity  0.00226;  drift  rate  +l.21°/day)  was  trimmed  to  a  nearly  synchronous  orbit  (semi¬ 
major  axis  42,159  km;  eeeentrieity  <  0.0001;  drift  rate  0.051  °/day).  The  thrust  obtained  from 
the  Gas  Thrusters  was  found  to  be  within  10  pereent  of  the  ground  measurements. 

The  Plasma  Thrusters  have  established  a  reeord  of  operation  in  orbit  for  eleetrie  propulsion 

systems.  The  measurement  of  their  performance  required  long  periods  of  operation,  and  thus 

elimination  of  the  natural  perturbations  from  the  motion.  This  is  beeause  the  effeet  of  the  Plasma 

Thrusters  is  quite  comparable  in  magnitude  to  that  of  the  geopotential  perturbation.  The  thrust 
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value  obtained  was  very  elose  to  the  value  measured  on  the  ground. 

The  quiet  period  from  12  November  to  22  December  was  neeessary  to  eheek  our  orbit  deter¬ 
mination  techniques  by  predieting  motion  under  the  effeet  of  the  natural  perturbations  only.  Our 
30-foot  antenna  and  receiver  installation  ean  measure  satellite  range  (via  time  delay  measure¬ 
ments  through  the  satellite  repeater)  very  aeeurately,  besides  measuring  azimuth  and  elevation. 
We  now  believe  that  our  orbit  determination  aeeuraey  is  of  the  order  of  0.3  km  in  semimajor  axis 
(typical  3a  value). 

The  Automatic  System  was  turned  on  for  a  long-term  experiment  on  2  3  December  1968.  At 
this  time  the  satellite  had  a  drift  rate  of  0.062  deg/day  west.  At  this  longitude  the  geopotential 
perturbation  tends  to  move  the  satellite  west.  The  Alarm  Command  was  sent  to  locate  the  station 
at  about  80°  West,  so  that  a  C-band  measurement  would  be  made  by  the  on-board  system.  This 
resulted  in  slowing  down  the  satellite,  as  seen  from  Fig.  22  and  recorded  by  the  orbit  determina¬ 
tions  of  17  January  and  20  February  (drift  rate  of  0.026  deg/day).  If  the  Automatic  System  had 
not  operated,  the  drift  rate  at  this  time  would  have  been  about  0.5  deg/day.  On  18  February  the 
station  was  moved  to  92.5°  West  longitude. 

The  data  obtained  from  the  sensors  were 
analyzed  to  test  several  hypotheses  about  their 
distribution.  Figure  2  3  shows  the  actual  un¬ 
smoothed  data  for  the  period  8  January  to  5 
March  1 969. 

To  test  the  validity  of  the  "weak-eoupling" 
approximation  to  the  drift  motion,  we  made 
least- square  fits  of  polynomials  of  time  to  the 
data.  The  results,  summarized  in  Table  11, 
show  that  the  unbiased  estimate  of  the  variance 
of  the  residuals  is  either  increased  or  not  sig¬ 
nificantly  decreased  in  going  from  a  second 
degree  to  a  third  degree  fit.  Thus  the  weak 
coupling  approximation,  whieh  predicts  a  near 
parabolic  behavior  for  small  longitude  excur¬ 
sions  (see  Eq.  (11-35)],  is  supported.  The  sec¬ 
ond  degree  least  square  fit  is  seen  in  Fig.  2  3 
(continuous  curve). 

Next  we  studied  the  residuals,  having  eliminated  the  trend  via  the  parabolie  least- square  fit. 

Stationarity:—  The  variance  of  the  residuals  is  different  in  the  different  periods,  as 
shown  in  Table  II.  Furthermore,  one  can  see  periods  in  which  the  data  scatter  is  greater,  even 
though  the  daily  thrust  has  not  changed  compared  to  a  relatively  quiet  period  (typically,  the 
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DAY  OF  OBSERVATION  (first  doy  equals  zero) 

Fig.  23.  Position  meosured  by  the  on-board  equipment 
in  a  typicol  period  of  operotion. 
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Fig.  24.  Autocorrelation  of  the  sensor 
errors  (first  period). 


Fig.  25.  Autocorrelation  of  the  sensor 
errors  (second  period). 
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TABLE  II 

VARIANCE  OF  RESIDUALS 

Linear 

Parabol  ic 

3rd  Degree 

4th  Degree 

Fit 

Fit 

Polynomial 

Polynomial 

Period 

(sec) 

(sec) 

(sec) 

(sec) 

First  (22  Oct. -12  Nov.) 

15.15 

11.42 

10.93 

11.59 

Second  (12  Nov. -15  Dec.) 

32.32 

28.87 

29.  16 

31.67 

Third  (23  Dec. -7  Jon.) 

61.98 

62.03 

63.05 

64.40 

Fourth  (8  Jan. -13  Feb.) 

39.63 

33.04 

32.26 

32.  16 

third  and  fourth  periods).  The  difference  in  variance  strongly  supports  the  hypothesis  that  the 
residuals  are  not  weakly  stationary  (and  thus  not  stationary)  in  time. 

Whiteness:—  We  applied  Anderson’s  modified  test  for  whiteness,  which  is  applicable 
if  the  residuals  are  uncorrelated.  This  is  done  by  computing  the  autocorrelation  for  unity  lag 
for  the  sample  from  the  time  series.  The  results  for  the  autocorrelation  function  are  shown  in 
Figs.  24  through  27.  We  also  applied  the  Durban  Watson  test  for  whiteness  to  the  different  periods. 
To  gain  some  intuitive  "feel"  for  the  behavior  of  the  autocorrelation  function  for  truncated  time 
series  obtained  from  a  "white”  process,  we  repeated  the  analysis  for  samples  obtained  from 


Fig.  26.  Autocorrelotion  of  the  sensor  errors 
(third  period). 


Fig.  27.  Autocorrelation  of  the  sensor  errors 
(fourth  period). 
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Fig.  28.  Autocorrelation  function 
for  white  Goussian  noise. 
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white  Gaussian  noise  and  uncorrelated  equiprobablc  noise.  A  typical  result  is  shown  in  Fig.  28. 
We  conclude  that  our  analysis  of  the  data  available  so  far  supports  the  hypothesis  that  the  resid¬ 
uals  are  a  white  random  process. 

llystograms:—  To  obtain  some  information  on  the  underlying  probability  function,  we 
prepared  hystograms  of  the  residuals.  The  hystograms  are  often  bimodal;  only  in  periods  111 
and  IV  do  the  hystograms  show  a  behavior  reminiscent  of  the  Gaussian  curve.  In  order  to  inter¬ 
pret  the  meaning  of  the  residuals,  it  is  important  to  note  that  the  telemetry  sampling  rate  is 
about  20  seconds.  Thus  there  is  an  intrinsic  uncertainty  of  about  0.1  degree  in  the  measurement. 
The  information  on  the  distribution  will  be  improved  as  the  available  measurement  sample  in¬ 
creases  in  time. 

VIII.  CONCLUSIONS 

The  Automatic  Orbit  Control  system  of  LAOS- 6  had  operated  successfully  for  several  months 
at  the  time  of  this  writing.  The  main  conclusions  to  be  drawn  from  this  experience  are,  in  our 
opinion: 

(1)  It  is  possible  to  automatically  control  synchronous  orbits  with  modest 
equipment  complexity.  This  leads  to  simplification  of  ground  operations 
and  reduction  in  operating  costs  of  networks  of  synchronous  satellites. 

(2)  The  accuracy  of  control,  which  in  the  present  experiments  appears  to  be 
of  a  few  degrees,  should  improve  in  future  systems. 

(3)  From  a  control- theoretic  point  of  view,  it  is  important  to  notice  the 
advantages  of  symmetric  measurements  and  thrusting,  which  ensure 
the  stability  of  the  coordinates  which  are  not  controllable  when  the  sat¬ 
ellite  longitude  is  sampled  once  per  day. 

The  data  presented  here  describe  the  first  few  months  of  operation  of  the  LES-6  automatic 
system,  using  the  Plasma  Thrusters.  Experiments  of  station  transfer  and  of  prolonged  opera¬ 
tion  with  the  Gas  Thrusters  will  be  carried  out  in  the  near  future. 


0.6 


WHITE  GAUSSIAN 
NOISE 


18-6-9786(1  ) 


*T  0.2 

< 

_l 
LU 

cr 
cr 
o 
o 

O  0 

=> 

<1 


Vi 


-0.2  -  O 


44 


In  the  author’s  opinion,  methods  of  Automatic  Orbit  Control  will  be  used  extensively  on 
tomorrow’s  synchronous  satellites.  This  judgment  is  based  on  the  faet  that  the  number  of  sat¬ 
ellites  in  synchronous  orbit  will  almost  eertainly  increase  greatly.  Stationkeeping  requirements 
will  then  be  dictated  by  the  need  to  avoid  interference  between  synchronous  satellites.  The  need 
for  Automatic  Orbit  Control  is  not  limited  to  the  synchronous  orbit.  Both  circular  nonsynchro- 
nous  and  general  elliptic  earth  orbits  are  perturbed  by  the  geopotential  harmonies,  the  sun  and 
the  moon.  The  system  that  has  been  flown  in  LES-6  is  only  a  beginning;  more  accurate  on-board 
control  of  synchronous  orbits  should  be  possible  on  the  next  generation  of  satellites.  We  also  see 
these  ideas  extended  and  applied  to  more  complex  problems  of  orbit  control. 
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APPENDIX  A 

THE  SOLAR  TIME  ON  AN  ARBITRARILY  ORIENTED  SATELLITE 


If  it  is  assumed  that  the  earth’s  and  satellite’s  spin  axes  are  parallel,  the  transit  of  the  sun 
on  corresponding  meridian  planes  of  the  two  bodies  occurs  at  times  inversely  proportional  to 
their  spin  rates.  This  is  no  longer  so  if  the  spin  axes  are  not  parallel.  The  variation  of  merid¬ 
ian  transit  time  has  two  major  causes:  the  inclination  of  the  plane  of  the  elliptic  to  the  equatorial 
body  of  the  earth  or  satellite,  and  the  elliptieity  of  the  earth's  orbit,  with  which  the  satellite’s 
orbit  around  the  sun  can  be  assumed  to  coincide  in  our  approximation. 

In  the  following,  formulas  are  given  for 
the  difference  between  mean  solar  time  and 
true  solar  time  measured  on  a  satellite  arbi¬ 
trarily  oriented  in  an  orbit  of  any  inclination, 
eccentricity,  and  period.  The  satellite  or¬ 
bital  inclination  does  not  affect  the  satellite 
solar  time  unless  errors  are  made  in  sensing 
the  earth.  The  formulas  take  the  earth's  orbit 
eccentricity  into  account  to  second  order,  and 
are  easily  extended  to  higher  order. 

The  reference  system  used  is  shown  in 
1*  ig.A-1.  The  unit  vectors  i,  j  are  in  the 
earth’ s  equatorial  plane,  and  i  coincides  with 
the  direction  of  the  earth’s  node.  The  triad 
i,  j,  k  is  right-handed.  The  ecliptic  plane  E  is 
defined  by  the  direction  cosines .  e{iij,  u^,  u^}; 
o  is  the  angle  between  the  ecliptic  plane  E 
and  the  satellite  equatorial  plane,  IT,  not  shown  in  the  figure;  ft  is  the  angle  between  plane  I 
and  the  orbital  plane,  e  and  e^,  are  the  eccentricities  of  the  satellite  and  earth  orbits,  respec¬ 
tively;  u:  27r  rad/sidereal  year;  fi  satellite  orbital  radian  frequency;  102.25°  argument 

of  perihelion  measured  from  the  earth's  nodal  line. 

The  coincidence  in  time  between  the  output  of  the  earth  and  sun  sensor  will  occur  whep  the 
projections  of  the  satellite  and  of  the  sun  differ  by  6,  which  is  the  angular  separation  between 
the  earth  and  sun  sensor,  and  an  integer  number  of  rotations. 

Then  the  times  of  successive  coincidences  are  found  to  be  the  solutions  of  the  equation 
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Fig.  A-l.  Satellite  reference  system. 
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No  loss  in  generality  occurs  if  one  takes  6  0,  which  displaces  the  time  reference. 

Let  us  now  take  as  a  reference  for  Eq.  (A-l )  that  at  which  the  first  coincidence  for  t  0, 
n  0  occurs  at  perihelion  for  a  satellite  spin  or  yaw  axis  normal  to  the  earth's  equatorial  plane. 
Then  u1  u2  0  and  u^  1  in  (A-2)  and  from  (A-l), 
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The  expression  (A-2)  is  substituted  into  (A-l).  The  times  of  coincidence  are  the  solutions 
of  Eq.  (A-l)  in  the  unknown  t,  for  different  integer  values  of  n. 

It  is  of  interest  to  consider  some  special  cases  and  to  rewrite  Eq.  (A-l)  in  ways  more  con¬ 
venient  for  the  applications. 

For  a  synchronous  satellite,  the  quantities  w  and  ft  can  be  expressed  in  rad/mean  solar 
day  as  w  -  27r|  ,  ft  =  2rr(l  +  4  ),  and  4  =  1/365.27.  Measuring  the  time  in  mean  solar  days,  it  is 
more  convenient  to  consider  the  angular  difference  A  between  the  satellite  and  sun  projection  on 
plane  U  at  the  end  of  each  successive  mean  solar  day.  Then, 

j  2e^  sin  (cot)  cos  (o) 

A  =  <p  —  tan  [cos  (o  )  tan  (2tt£  n  +  <p,  )]  -  - - - ^ ~ - -  ,  (A -3) 

1  —  sin  (o' )  sin  (cot) 

where  <p  is  given  by  the  right  side  of  (A-l)  with  ft  =  2 7r ( 1  +  4  )  and  with  <P SQ  given  by  (A-2). 

For  a  synchronous  equatorial  satellite  in  a  circular  orbit  with  its  spin  or  yaw  axis  perpen¬ 
dicular  to  the  equatorial  plane  of  the  earth,  Eq.  (A - 3 )  yields  an  approximate  expression  for  the 
sun  transit  time  correction  on  the  earth  Ap,  being 

p  =  0  4  ~  1/365.27 

f  =  0  ar  102.25° 

U1  =  u2  '  0  ’  fE  “  0.01675 

o  =  tj  =  23.45°  ,  u^  =  1 

2<.  sin  (2tt4  n)  cos  (77  ) 

A,  27t|  n  +  if  -  tan  [cos  (77  )  tan  ( 2tt ^  n  +  </?,)]  -  - 1 — ^ ^ - .  (A -4) 

so  1  -  sin4  (77  )  sin4  (2tt|  n) 

For  a  synchronous  equatorial  satellite  in  a  circular  orbit  with  its  spin  or  yaw  axis  tilted 
around  the  earth’s  line  of  nodes,  the  foregoing  values  are  used,  except  that,  4  being  the  angle 
of  tilt, 


u2  sin  {4  ) 


cos  {a )  sin  (77  )  sin  (4)  +  cos  (7 7  )  cos  {4  )  cos  (77  -4) 


and 


2f  sin  (27 1-4  n)  cos  (77  -  4  ) 

A 4  2tt|  n  +  ip  -  tan  [cos  (77  -4)  tan  (27r|  n  +  <p,  )]  -  - - - x - ^ —  -  .  (A -5) 

1  -  sin4  (77  —  4)  sin4  (27:4  n) 


Comparing  this  equation  with  (A -4).  one  sees  that  a  rotation  of  the  yaw  axis  around  the 
earth’ s  line  of  nodes  is  equivalent  to  a  rotation  of  the  plane  of  the  ecliptic  around  the  same  line. 
This  means  that  the  satellite  solstices  and  equinoxes  take  place  on  the  same  day  as  the  earth’s, 
as  is  obvious  from  the  fact  that  the  lines  of  nodes  of  satellite  and  earth  still  coincide.  Plots  of 
angular  displacement  vs  time  of  the  year  for  this  case  are  shown  in  Fig.  A-2. 

For  a  synchronous  equatorial  satellite  in  a  circular  orbit  with  its  yaw  axis  tilted  around  the 
normal  to  the  earth’s  line  of  nodes,  and  with  x  being  the  angle  of  tilt, 

u1  sin(x)  .  u2  0  .  u^  cos(x)  .  cos  (cv )  cos  (x )  cos  (77  ) 
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Fig.  A-2.  Salar  time  as  seen  an  the  sotellite 
for  different  angles  of  tilt. 


and  therefore 

A  =  +  2?r£  n  4  <p  —  tan  [eos  (77  )  cos  (x)  tan  (2tt£  n  +  </>,  )| 

X  so  n 

2e^  sin  (2tt£  n)  cos  (x)  eos  (77  ) 

2  2  2  2  *  (A_h) 

eos  (2tt£  n)  +  eos  (77  )  eos  (x)  sin  ( 2tt ^  n) 

From  (A-l)  and  U  2n(i  +  £  ),  the  eccentricity  of  the  satellite  orbit  is  seen  to  contribute 
the  term 

Ar  2c  sin[27r(l  +  $  )  n]  eos  (ft)  /A 

^  - * - —x -  .  [J\  -  7  j 

1  sm  (0)  sin  [2t(1  +  {  )  n] 

This  confirms  that  the  eccentricity  effect  can  be  eliminated  to  second  order  in  c  by  taking 
two  measurements  of  satellite  solar  time  at  times  n  and  n  +  (1/2),  being 

A'(n)  ~  -  A'[n  +  (i)]  , 

since  £  is  small. 
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APPENDIX  B 

THE  OBSERVABILITY  AND  CONTROLLABILITY  MATRICES 
FOR  THE  SAMPLED  CASE  + 

The  sampled  case  equations  are  given  by 

x(kT  +  T)  -  A^xfkT)  +  Bgu(kT)  ,  (B-l) 

and  the  measured  vector  is 

y(k)  Cx(k)  ,  (B-2) 

where  the  matrices  and  are  given  by  Eqs.  (IV-20),  (IV-21),  (IV-22)  for  the  different 
cases. 

The  system  is  again  factored  into  two  parts  and  we  can  examine  both  separately. 

i.  Case  T  =  27r/co  (sampling  once  per  orbit) 

Ifu2  0, 


0  0  0  0 

1111 

0  0  0  0 

0  0  0  0 


(H-3) 


This  is  singular  so  that  the  sampled  system  is  not  controllable  from  u^.  (Actually,  we  could 
have  inferred  this  from  the  continuous  case.) 

Let  us  check  the  case  T  27r/co  and  u^  0.  Here  we  have 


~2 
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-  1  87T 
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-  247T 


(B-4) 


and  we  see  that  sampling  at  the  rate  T  27r/co  has  destroyed  controllability  from  u^. 

If  we  wish  to  check  controllability  from  both  u^  and  u^  together,  the  relevant  matrix  is 


t  From  Ref.  27. 
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0  —  1  2tt  0 
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CO 


•  24tt 
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which  is  also  singular.  Hence,  even  with  both  and  available  one  cannot  control  if  only 
one  thrust  per  orbit  is  used. 

The  controllability  of  x^  and  x ^  from  u^  is  revealed  by 


!<4 


0  0 

1  1 


(H-6) 


which  is  singular  and  hence  indicates  a  lack  of  controllability. 

The  observability  picture  is  as  follows.  From  x^  we  have  an  observability  matrix 
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From  x^  we  have  an  observability  matrix  of  the  form 


—2 


0  0  1 


-  1  2  7T  0  1 


24tt  0  1 


—  36tt  0  1 


—  67 r 


1  Zn 


18tt 


(B-7) 


(B-8) 


and  from  x^  and  together  we  have 
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—  1 27 r  0  1  — — 
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Clearly  none  of  these  matrices  is  of  rank  4  so  we  cannot  observe  the  inplane  components  of  the 
motion  using  one  observation  per  orbit. 

The  out-of-plane  situation  is  revealed  by 


—4 


1  0 
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so  we  see  that  this  motion  is  not  observable  either. 
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ii.  Case  T  tt/co  (sampling  twice  per  orbit) 

In  this  case,  the  controllability  matrices  K  ^ and  take  the  form 
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Both  and  are  singular. 

but  K3 

is  of 

rank  4  (the  firs 

t  four 

columns  are 

independent).  Hence 

it  is  possible  to  control  the  in 

-plane 

motion  by  thrusting  only  twice  per 

orbit 

,  provided  that  one 

has  both  tangential  and  radial  thrusting  capability. 

For  the  out-of-plane  motion,  the  situation  is  not  improved,  however,  because 
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and  is  not  of  rank  2. 

As  far  as  observability  is  concerned,  we  have 
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Only  matrix  is  of  rank  4  (the  first  four  rows  are  independent).  Hence  the  in-plane  motion 
is  observable,  sampling  only  twice  per  orbit  if  both  angle  and  altitude  ran  be  measured. 

1  or  the  out- of- plane  motion,  we  need  to  examine 
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and  we  see  that  this  part  of  the  system  is  still  not  observable. 

iii.  Case  T  rj  2oj  (sampling  four  times  per  orbit) 

The  controllability  matrices  K^,  K  },  and  are  now  given  by 
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K  and  K ^  arc 

?  of 

rank  4  but  K 

.  is  not. 

1 

Thus,  controlling  the 

()  dependent 

s  one  control  the'  e 

nt ire  state*  if 

thrusting  is  done*  four  time 

“S  per  orbit. 

For  the  out-of-plano  motion,  takes  the  form 
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1  bus  this  motion  is  controllable. 

To  check  observability  in  the  four  times  per  orbit  sampled  case,  we  need 
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Clearly,  is  not  of  rank  4  so  we  cannot  observe  the  state  by  looking  at  altitude  alone.  How¬ 
ever,  I_  and  L ^  are  of  rank  4,  so  observation  from  O,  and  from  r  and  0  together,  is  possible. 
The  out- of- plane  motion  is  covered  by 
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and  it  is  observable  from  four  samples  of  cp  per  orbit. 
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